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This  report  was  originally  prepared  under  the  sponsorship  of  the 
U  S.  Navy  through  the  Office  of  Naval  Research  (Mechanics  Branch), 
Contract  No.  2653(00),  and  was  presented  at  the  IAS  31st  Annual  Meeting, 
New  York,  New  York,  January  21-23,  1963. 

As  pointed  out  in  the  text,  a  part  of  the  study  with  the  wedge  and  cone 
problems  was  carried  out  for  the  U.  S.  Air  Force  through  the  Office  of 
Scientific  Research,  under  Contract  AF  49(638)-952. 

Since  the  completion  of  the  original  version  of  the  report,  which 
has  been  distributed  as  IAS  paper  No.  63-92,  errata  as  well  as  points  of 
ambiguity  have  been  called  to  the  writer's  attention;  this  version  has  been 
revised  to  incorporate  the  desirable  changes. 
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THE  BLUNT  -BODY  PROBLEM  iN  HYPERSONIC  FLOW 
AT  LOW  REYNOLDS  NUMBER 

H  K  CHENG 

Cornell  Aeronautical  Laborator/  Inc,,  Buffalo,  New  York 

ABSTRACT 

Existing  theoretical  analyses  and  experimental  results  of  the  stagnation 
region  in  hypersonic  flow  at  low  Reynolds  number  are  discussed  The  approach 
based  on  the  thin- shoe k -lay er  appro imaticn  is  extended  to  study  flow  fields 
beyond  the  stagnation  region  as  well  as  in  the  shock- transition  zone. 

The  basic  flew  model  consists  cf  two  adjoining  thin  layers’  a  shock- 
transition  zone  and  a  shock  layer.  The  sys;  m  of  partial  differential  equations 
governing  the  high-density  shock  layer  reduces  to  the  parabolic  type.  The 
system  governing  the  shock  •  transition  zone  reduces  to  ordinary  differential 
equations  similar  tc  those  cf  the  one  -  dimensional  shock  wave.  They  give  rise 
to  a  set  of  conservation  relations  across  the  shock,  which  account  for  the  trans¬ 
port  processes  immediately  behind  the  shock  but  do  not  involve  details  of  shock- 
wave  structure.  With  the  modified  Rankine-Hugoniot  relations,  the  flow  field 
in  the  shock  layer  can  be  determined  independently  of  the  shock-transition  zone. 
An  essential  feature  of  this  formulation  is  that  when  applied  in  conjunction  with 
nonslip  surface  conditions  it  always  yields  the  appropriate  surface  heat-trans¬ 
fer  rate  and  skin  friction  (for  unit  acc  ommedation  coefficients)  in  the  free- 
molecule  limit  Com  pan  son  with  more  exact  solutions  of  Levinsky  and 

Yoshihara,  and  of  Van  Dyke  reveals  the  adequacy  ot  the  present  theory  in  both 
the  merged-layer  regime  and  the  mviscid  limit 

A  finite-difference  method  is  developed  for  solving  the  shock-layer 
equations.  As  examples  of  llcw  around  simple  blunt  or  nonslender  shapes,  the 
problems  of  the  paraboloid  hyperboloid  as  well  as  the  wedge  and  cone  have 
been  analyzed  for  a  perfect  gas  with  constant  specific  heats.  Generally,  solu¬ 
tions  for  low  Reynolds  numoer  reveal  substantial  changes  in  total  enthalpy  and 
tangential  velocity  across  the  shock,  analogous  to  surface  slip  effects.  The 
solutions  provide  smooth  transitions  from  the  boundary-layer  theory  to  the  free- 
molecule  limit.  The  velocity  and  temperature  gradients  near  the  inner  edge  of 
the  shock  layer  increase  with  distance  from  the  stagnation  point  or  the  apex, 
indicative  of  boundary-layer  development.  Results  obtained  for  paraboloids 
and  hyperboloids  indicate  the  existence  cf  a  remarkably  wide  range  of  Reynolds 
number  (corresponding  to  the  boundary-layer  through  the  mcipient-merged- 
layer  regime',  in  which  the  distributions  of  heat-transfer  rate  and  skin  friction 
on  a  smooth  blunt  body  remain  essentially  unchanged. 

Also  studied  is  nonequilibrium  dissociation  with  species  diffusion  in 
shock  layers  under  influence  of  strong  surface  cooling.  Solutions  for  the 
stagnation  region  based  on  -a  single  dis sociation- recombination  kinetic  model 
show  a  substantial  concentration  jump  across  the  shock  and  extensive  influ¬ 
ence  of  surface  cooling  ai  low  Reynolds  number  The  role  of  pressure  diffusion 
in  a  viscous  shock  layer  is  also  discussed 
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pertaining  to  the  stagnation  condition  behind  the  bow  shock  in 
an  inviscid  flow 
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INTRODUCTION 


A  number  of  effective  methods  have  been  developed  for  analyzing 

inviscid  hypersonic  flows  and  the  associated  boundary  layers  around  blunt 

bodies  (see,  for  example.  Ref.  1-4).  Much  less  has  been  done,  however, 

for  the  corresponding  problems  at  low  Reynolds  number  or  low  density. 

The  analysis  of  the  low-density  phenomena  is  nevertheless  basic  to  the  study 

of  many  problems  of  high-speed  flight  at  high  altitude  where  the  thin-boundary- 

layer  concept  no  longer  applies.  In  this  paper,  the  problem  of  hypersonic 

rarefied- gas  flow  will  be  studied  within  the  framework  of  the  conlinuum  theory. 

With  an  increasing  degree  of  rarefaction,  the  continuum  theory  begins 

to  break  down.  Consequently,  a  valid  description  of  the  flow  field  will  have 

to  be  based  on  the  kinetic  theory  to  account  for  the  collision  processes.  The 

transition  from  the  continuum  to  the  free-molecuie  limits  from  the  view- 

5 

point  of  flight  application  was  considered  some  time  ago  by  Tsien.  Tsien 
classified  the  problems  of  rarefied-gas  flow  into  four  regimes:  the  continuum, 
the  slip-flow,  the  transitional,  and  the  free -molecule  regimes.  It  was  pro¬ 
posed  at  that  time  that  departure  from  the  classical  continuum  model  in  the 
slip-flow  regime  may  be  accounted  for  by  including  in  the  governing  equations 
the  higher-order  (Burnett)  terms  in  a  formal  development  for  the  solution  to 
the  Boltzmann  equation.  ^  However,  the  growing  body  of  evidence  from  both 
experimental  and  theoretical  studies  has  since  shown  disagreement  with 
these  higher-order  equations,  and,  in  most  instances  has  indicated  that  the 
relatively  simpler  Navie  r -Stroke  s  equations,  with  appropriate  boundary 

conditions  to  account  for  the  slip  effects  (and  transport  coefficients  evaluated 

7  8 

according  to  the  Chapman-Enskog  theory),  are  superior.  ’  A  continuum- 
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flow  model  based  on  the  Navier-Stokes  equations  thus  promises  a  domain 

of  applicability  far  more  extensive  than  previously  conceived  This 
9-12 

view  is  shared  by  the  writer,  and  the  present  study  attempts  to  extend 

the  continuum-flow  theory,  as  far  as  is  permissible,  into  the  range  of  rarefied 
gas  dynamics. 

It  is  apparent  that  solution  tc  the  transitional  regime,  as  well  as  the 

final  justification  of  the  continuum  analysis,  will  have  to  come  from  the 

Boltzmann  equation  governing  the  velccity-distribution  function.  A  number 

of  approximate  methods  have  been  developed  in  recent  years  to  overcome 

13  14  1  5 

the  difficulty  of  solving  the  Boltzmann  equation.  '  ’  However,  a  method 

which  is  satisfactory  in  both  accuracy  and  simplicity  has  yet  to  beccme 
available  Of  significance  in  this  respect  is  the  recent  study  of  the  shock- 
transition  zone  by  Liepman  and  co-workers,  in  which  a  variant  of  Krook's 
collision  model*  ^  is  adopted.  The  solution  obtained  shows  little  difference 
from  the  Navier-Stokes  solution  on  the  higher-pressure  side  of  the  shock 
even  for  a  Mach  number  of  ten.  It  seams  to  provide  support  to  the  idea  that 
the  Navier-Stokes  equations  are  adequate  for  describing  flow  fields  even  when 
the  shockwave  structure  becomes  important. 

The  problem  of  hypersonic  flow  cf  rarefied  gas  around  a  blunt  body 
has  been  discussed  previously  by  Probstein  in  Ref.  11  where  an  extensive 


With  this  model,  the  Boltzmann  equation  is  simplified  to  one  formally 
resembling  a  relaxation  equation 

To  further  augment  this  idea,  it  would  be  desirable  to  include  in  the 
shock- structure  analysis  of  Liepman  and  co-workers  the  effects  of  shock 
curvature  as  well  as  velocity  and  temperature  gradients  behind  the  shock. 
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account  of  the  analyses  of  the  stagnation  region  in  the  various  flow  regimes 

1  2 

is  also  given.  Recently,  the  problem  has  been  reviewed  by  Van  Dyke, 
err  pna  si  zing  the  second-order  boundary-layer  theory.  The  present  report 
consists  of  two  parts:  In  Part  I,  the  existing  analyses  of  the  stagnation 
region  will  be  discussed  with  emphasis  on  the  more  recent  studies  related 
to  the  viscous-layer  and  merged-layer  regimes.  A  comparison  with 
experimental  results  is  included.  The  discussion  is  oriented  toward  the 
thin- shock- layer  approach  which  is  to  be  further  developed  in  Part  II. 

In  Part  II,  the  thin- shock-layer  approach  to  the  viscous  blunt-body 
problem  is  extended  to  study  flow  fields  other  than,  or  removed  from, 
the  stagnation  region.  The  new  development  includes  an  extension  of  the 
thin-layer  approximation  to  analyze  the  shock-transition  zone  and  a  finite  - 
difference  method  for  solving  the  shock  layer  equations  without  similarity 
assumptions.  Among  the  particular  problems  studied  are  viscous  hyper¬ 
sonic  flows  over  wedges,  paraboloids,  and  other  simple  nonslender  shapes, 
as  well  as  nonequilibrium  dissociation  with  species  diffusion  in  the  shock 
layer  under  the  influence  of  strong  surface  cooling. 
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PART  I.  REVIEW  OF  EXISTING  RESEARCH  ON  STAGNATION  REGION 


The  manner  in  which  flow  regimes  are  classified  is,  of  course,  quite 
arbitrary,  and  varies  from  one  problem  to  another.  For  a  stagnation  region 
in  hypersonic  flow,  Probstein  and  Kemp  observe  that  there  are  six  regimes 
between  the  boundary-layer  and  free-mclecule  limits.  Figure  1,  reproduced 
from  Ref.  11,  illustrates  the  extent  of  the  various  regimes  in  terms  of  flight 
speed  and  altitude  The  boundaries  delimiting  each  of  these  regimes  were 
originally  estimated  by  Probstein^  ^  for  an  ARDC  1959  model  atmosphere 
and  a  cooled  body  of  one -foot  nose  radius,  assuming  the  compression  ratio 
across  the  shock  to  be  ten.  These  regimes  numbered  below  in  the  order  of 
increasing  degree  of  rarefaction,  may  be  characterized  by  the  following 
de  scriptions . 

(1)  Vorlic ity-Interaction  Regime  -  The  boundary  layer  is  still  thin, 
but  the  second-order  effects  become  appreciable.  They  include  slip  and 
temperature  jump,  surface  curvatures  and  the  interaction  of  the  boundary 
layer,  through  its  finite  thickness,  with  the  outer  flow.  For  the  axisym- 
metric  case,  the  most  important  of  these  effects  is  that  resulting  from  the 
external  vorticity. 

(2)  Viscous-Layer  Regime  -  The  viscous  region  becomes  an  appreciable 
fraction  of  the  shock  layer  and  may  be  too  thick  to  permit  application  of  the 
boundary-layer  concept. 

(3)  Incipient-Merged  Layer  Regime  -  The  viscous  and  other  transport 
effects  are  no  longer  negligible  behind  the  shock,  and  the  thickness  of  the  shock 
itself  becomes  an  appreciable  fraction  of  the  shock  layer.  Hence  the  shock  as 
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a  surface  of  discontinuity  may  no  longer  be  an  adequate  description. 

(4)  Fully-Merged  Layer  Regime  -  The  shock  thickness  is  comparable 
to,  and  perhaps  even  greater  than,  that  of  the  shock  layer.  Not  only  is  the 
shock  layer  fully  viscous,  but  the  shock-transition  zone  may  no  longer  be  dis¬ 
tinguished  from  other  parts  of  the  flow  field. 

(5)  Transitional  Regime  -  The  collision  frequency  among  gas  particles 
is  not  high  enough  to  make  the  continuum  description  valid,  yet  not  sufficiently 
low  to  make  the  kinetic-theory  description  simple. 

(6)  First  Collision  Regime  -  There  is  only  a  small  fraction  of  incident 
particles  experiencing  one  collision  or  more  before  reaching  the  body 
surface . 

The  pre  sent  report  mainly  concerns  regimes  (l)-(4).  The  study  will  be 
made  under  the  assumption  of  a  continuum  flow  model  governed  by  the  Navier- 
Stokes  equations.  That  assumption  has  in  fact  been  implicit  in  the  determination 
of  the  various  boundaries  in  Fig.  1  .  For  later  reference,  the  regimes  (3)  and 
(4)  together  may  be  called  the  merged-layer  regime. 

Aside  from  providing  a  general  picture  of  the  flow  field  at  various  degrees 
of  rarefaction,  the  above  classification  serves  to  indicate  the  appropriate  flow 
models  for  analysis.  Attempts  to  analyze  the  stagnation  region  for  regimes 
(l)-(4)  have  resulted  in  three  main  approaches.  One,  which  applies  to  regime 
(1),  assumes  small  departure  from  the  boundary- layer  limit.  This  approach, 
which  has  been  most  extensively  treated,  may  be  identified  with  the  higher- 
order  boundary-layer  theories.  In  another,  applicable  to  regimes  (2),  (3), 
and  partly  to  (4),  one  considers  viscous  and  other  transport  effects  throughout 
the  shock  layer,  and  integrates  numerically  a  system  of  simplified  Navier- 
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Stokes  equations  along  the  axis  of  symmetry.  Encompassing  both  of  these 

approaches  is  a  third  approach,  favored  by  the  present  author,  -which  utilizes 

the  thin- shock-layer  approximation  familiar  from  the  invi  scid  hyper  sonic 
18-21 

flow  theory.  The  analysis  based  on  this  third  approach  oniy  requires 

the  distinction  between  two  flow  regimes,  labelled  I  and  II  in  Fig.  1.1.  Regime 
I  contains  regimes  {])•  and  (2)  of  the  foregoing  classification;  Regime  II  contains 
(3)  and  the  main  part  of  (4).  Current  theories  will  be  discussed  under  these 
three  main  approaches  in  Sections  t-3.  The  results  of  theories  and  experiments 
from  several  sources  are  compared  in  Section  4. 

1.  SECOND-ORDER  BOUNDARY  - LA  V ER  THEORIES 

The  second-order  corrections  to  the  boundary-layer  theory  may  be  sub¬ 
divided  in  various  categories:  the  external  vorticity  effect,  slip  and  temperature- 
jump  effects;  the  non-vanishing  suriace  curvature  effects;  and  the  boundary-layer 
displacement  effect.  The  boundary-layer  displacement  affects  directly  the 
matching  of  the  boundary  layer  with  the  outer  vortical  flow,  and  also  induces 
pressure  and  velocity  changes  in  the  outer  ''mviscid  flow"  On  account  of 
their  subtle  nature,  the  subdivision  concerning  displacement  effects  has  been 
rather  arbitrary. 

4 

The  relatively  earlier  analyses  of  Hayes  and  Probstein,  Rott  and 

Lenard^  and  Kemp^  are  devoted  primarily  to  the  effect  of  external  vorticity. 

More  systematic  treatments  ot  the  second- order  effects  are  subsequently  given 
24  25  2b 

by  Van  Dyke,  Lenard  and  Maslen  The  analysis  of  Ferri.  Zakkay  and 
27 

Ting  considers  vorticity  effect  and  may  also  be  regarded  as  belonging  to  this 
class.  There  has  been  considerable  disagreement  among  the  various  analyses 
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of  slip  and  temperature -jump  effects,  as  well  as  the  vorticity  and  displacement 

1  2 

effects.  Readers  are  referred.to  Van  Dyke's  review  for  the  many  critical 
discussions  and  comments.  The  following  discussion  will  concern  only  certain 
outstanding  aspects  of  the  higher-order  theories  which  must  be  considered  sub¬ 
sequently  in  this  paper:  the  external  vorticity  effect,  the  displacement  effect, 
and  the  slip  and  the  temperature-jump  effects. 

1.1  External  Vorticity  Effect 

28 

In  his  earlier  work  on  the  external-vorticity  effect,  Hayes  assumes 
that  the  boundary-layer  displacement  dees  not  appreciably  change  the  pressure 
and  velocity  distributions  along  each  streamline  in  the  outer  flow.  He  arrives 
at  a  modified  outer  boundary  condition  for  the  velocity  field  f  in  the  form 

=  /  +  2  Pi  f  1 - >  o<?  (1.1) 

where  I~L  ,  termed  the  vorticity -interac tion  parameter,  is  the  ratio  of  the 

external  vorticity  to  some  average  value  of  vorticity  in  the  boundary  layer. 

This  condition  was  applied  to  study  the  axisymmetric  stagnation-point  boundary 

4 

layer  by  Hayes  and  Probstein.  In  this  case,  T2-  is  given  by  Hayes  and 
Probstein  as 


XL  = 


(p  /p  V,_!i L) 

V  wyrz  )y  3  / 


(1.2) 


yr  e  /nr  '  p*0**- 

Equation  (1.1)  shows  clearly  that  external  vorticity  effectively  increases  the 
velocity  at  the  boundary  layer's  outer  edge.  One  may  note  that  since  Eq.  (1.1) 
does  not  contain  r|  ,  it  also  allows  for  the  streamline  displacement  in  the 
matching  of  the  inner  and  outer  velocity  fields. 

The  above  boundary  condition  may  be  replaced  simply  by 

(1.3) 


rl*f 
d  nz 


->  a 


ft 


->  oo. 
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The  equivalence  of  Eqs.  (1  1)  and  (i  .  3)  is  guaranteed  by  the  differential 

23  29 

equations  governing  f  >  as  n°te<T  by  Kemp  and  Moore 

The  method  based  on  the  local-similarity  principle  has  been  a  very 

useful  tool  for  analysis  of  heat  transfer  from  a  hypersonic  boundary  layer  to 

a  relatively  cold  bedy.^’  ^  It  is  natural  to  anticipate  an  extension  of  the 

method  to  account  for  the  external- vorticity  effect.  This  is  presumably  the 

27 

viewpoint  of  Ferri,  Zakkay  and  Ting  who  treat  the  vorticity  interaction  on 
the  basis  of  a  local-flat-plate  similarity.  This  basis  implies  however  that 
the  velocity  profile  is  completely  unaffected  by  the  tangential  pressure  gradient 
as  well  as  by  the  e.xternal  vorticity.  The  outer  inviscid  flow  is  assumed  to  have 
a  linear  velocity  profile  and  is  completely  unaltered  (and  undeflected)  by  the 
presence  of  the  boundary  layer.  The  procedure  of  joining  the  inner  and  outer 
solutions  appears  to  be  a  rather  complicated  one  which  involves  matching  the 
velocity  and  temperature  fields  ever  two  separate  boundaries.  One  of  the 
final  equations  used  (Eqs.  (341  or  (38)  of  Ref.  27)  may  nevertheless  be  iden¬ 
tified  in  form  with  Hayes'  boundary  condition  Eq.  (1  .  1). 

1,2  Displacement  Effect 

Except  for  the  displacement  effect  implicit  in  the  boundary  condition 

of  Eq.  (1.1),  Hayes  and  Probstein's  vc  rticity-mteraction  theory  assumes  no 

displacement  effect  on  the  outer  flow,  This  assumption  is  justified  in  the 

case  of  a  thin  shock  layer.  In  order  to  account  for  the  displacement  effect 

completely,  one  has  to  solve  the  difficult  elliptic  problem  (presumably  by 

1  2 

iteration),  and  this  has  been  done  only  by  Van  Dyke.  Unfortunately,  the 

problem  is  further  obscured  by  an  apparent  confusion  over  the  division  of 

24 

the  vorticity  and  displacement  effects.  While  Van  Dyke's  analysis  is 
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believed  to  be  the  most  systematic  and  elegant,  his  subdivision  of  the  inter¬ 
action  problem  into  the  effect  of  "entropy  gradient"  (i.e.  vorticity)  and  the 
effect  of  "displacement  speed"  is  not  completely  satisfactory  In  each  of 
these  effects,  according  to  this  classification,  there  is  a  second-order  effect 
appearing  in  both  the  outer  boundary  condition  as  well  as  in  the  governing 
differential  equations.  The  second-order  corrections  in  the  differential 
equations  actually  arise  from  the  change  in  the  pressure  at  the  boundary 
layer's  outer  edge,  and  may  therefore  be  combined  as  the  displacement 
effect  on  the  external  pressure.  The  combined  effect  turns  out  to  be,  in 
most  cases,  much  smaller  than  if  separately  calculated,  and,  as  already 
noted,  is  negligible  in  the  thin- shock- layer  approximation. 

1  .  3  Slip  and  T emperature-Jump  Effects 

22  25  24 

While  the  analyses  of  Rott  and  Lenard,  ’  as  well  as  Van  Dyke, 
show  sizable  slip  effects,  Maslen's  analysis^  reveals  that  these  effects 
are  definitely  negligible  for  heat-transfer  analysis.  It  may  be  noted  that  the 
results  of  these  analyses  depend,  of  course,  on  the  specific  coefficients 
entering  the  slip  and  temperature-jump  conditions.  The  large  temperature- 
jump  effects  reported  by  Rott  and  Lenard,  and  by  Van  Dyke  are  rather  sur¬ 
prising,  since  these  effects  are  known  to  vanish  in  either  the  limit  of  a 

vanishing  wall-to- stagnation  temperature  ratio  or  the  limit  of  a  vanishingly 
1  2 

thin  shock  layer. 


The  change  in  velocity  at  the  outer  edge  associated  with  the  boundary- 
layer  displacement. 
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2.  ANALYSIS  OF  VISCOUS  SHOCK  LAYERBASED  ON  SIMILARIT  Y  ASSUMPTION 


When  departure  from  the  boundary-layer  limit  is  too  large  for  the  con¬ 
cept  to  be  applicable,  the  full  Navier -Stoke s  or  other  appropriate  equations 
must  be  used  to  describe  the  stagnation  region.  Most  of  the  existing  analyses 
in  this  class  use  the  Rankine -Hugoniot  shock  relations  as  an  outer  boundary 
condition,  and,  therefore,  can  be  applied  only  to  the  viscous-layer  regime. 
However,  in  the  merged-layer  regimes,  the  Rankine-Hugoniot  relations  are 
not  applicable  and  have  to  be  modified.  In  order  to  make  the  analyses  tract¬ 
able,  it  has  been  found  necessary  to  invoke  the  similarity  assumption  for  the 
flow  field  in  the  stagnation  region.  This  assumption  is  justified  for  two  types 
of  models:  one  assumes  a  constant  density,  the  other  a  thin  shock  layer, 

2,  1  Similar  Solutions  Based  on  Constant  Density 

Assuming  a  constant  density  behind  a  spherical  cr  cylindrical  shock  which 

l  o 

is  concentric  with  a  spherical  or  cylindrical  body,  Probstein  and  Kemp'”  have 
reduced  the  Navier -Stoke s  equations  in  the  stagnation  region  to  a  system  of 
ordinary  differential  equations.  Solutions  are  obtained  for  a  cold  surface  both 
in  the  vi scous- layer  and  in  the  incipient-merged-layer  regimes.  The  analysis 
in  the  latter  regime  must  include  the  thicKness-curvature  effect  of  the  shock 
and  the  modification  of  the  Rankine-Hugoniot  relations  to  allow  for  the  trans¬ 
port  effects  immediately  behind  the  shock.  In  Probstein  and  Kemp's  model,  ^ 
the  shock-transition  zone  is  treated  separately  from,  but  simultaneously 
with,  the  shock  layer.  In  spite  of  the  many  idealizations  of  the  flow 

The  similarity  assumption  amounts  to  postulating  the  flow  field  to  be  in  a 
certain  form,  so  that  a  separation  of  the  variables  may  be  achieved  in  the 
governing  differential  equations.  By  the  similarity,  the  problem  is  reduced 
to  one  involving  only  ordinary  differential  equations. 
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model,  the  analysis  of  Probstein  and  Kemp  was  first  to  reveal  the  many 
important  features  of  the  viscous-layer  and  incipient-mer ged-layer  regimes. 
This  constant-density  model  perhaps  describes  more  appropriately  the 
physical  problem  of  an  insulated  body,  which  has  been  considered  by  Hoshizaki" 

in  the  vi scous -layer  regime.  The  constant- density  viscous-layer  model  has 

33  34 

been  analyzed  by  Oguchi  to  include  a  magnetic  field.  Oguchi  has,  in  fact, 


applied  a  systematic  expansion  in  the  density  ratio 


to  dc  scribe 


analytically  the  constant-density  vi  scous -layer  model.  Although  the  range  of 

validity  of  Oguchi's  solution  is  limited  by  the  viscous -layer  model  and  by  the 

breakdown  of  the  formal  expansion  at  high  Reynolds  number,  the  simplicity  to 

be  gained  by  the  thin- shock-layer  approach  is  quite  evident  from  the  analysis. 

2,  2  Similar  Solution  Based  on  Thin  Shock  Layer 

The  stagnation-point  boundary  layer  is  much  easier  to  handle  than  the 

corresponding  problem  based  on  the  full  Navier-Stokes  equations  because  the 

boundary-layer  equations  are  parabolic  and  amenable  to  separation  of  variables 

for  the  stagnation-point  problem.  The  Navier-Stokes  equations  may,  however, 

be  simplified  also  to  the  parabolic  type  by  assuming  that  the  shock  layer  (which 

includes  the  boundary  layer  as  a  part)  is  thin.  In  this  way,  the  similarity 

assumption  in  the  stagnation  region  can  be  justified  without  the  assumption  of 

34 

constant  density.  Thus,  Ho  and  F  robstein  integrate  numerically  the 
simplified  Navier-Stokes  equations  along  the  axis  of  symmetry  from  the  body 
to  the  shock  in  the  viscous-layer  regime.  The  reduced  equations  are  similar 
to  the  boundary-layer  equations  except  that  the  outer  boundary  is  now  at  a 
finite  distance  from  the  surface  (where  the  Rankine-Hugoniot  shock  conditions 
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are  to  be  satisfied)  and  that  the  tangential  pressure  gradient  is  not  a  constant, 

35 

Herring  postulated  simply  that  the  boundary-layer  equations  hold  between 
the  body  and  the  shock,  His  analysis  of  the  stagnation  region  is  then  virtually 
the  same  as  Ho  and  Probstein's,  except  for  the  fact  that  the  tangential  pres¬ 
sure  gradient  term  in  the  momentum  equation  is  a  constant  determined  by  its 

36  v 

value  at  the  shock.  Chung  studied  nonequilibrium  dissociation  in  the  viscous 

stagnation  region.  He  carried  out  his  study  within  the  framework  of  the  viscous- 

layer  model  of  Ho  and  Probstein,  but  evaluates  the  tangential-pressure- 

gradient  term  as  in  Herring's  analysis.  The  results  of  Chung's  calculation 

will  be  discussed  later  in  connection  with  the  study  of  Section  7  of  Part  II. 

37 

Also  using  a  viscous-layer  model.  Goldberg  and  Scala  studied  mass  trans¬ 
fer  in  the  hypersonic  shock  layer.  The  validity  of  some  of  their  results  may 
be  questionable  because  of  the  local  chemical  equilibrium  postulated. 

The  use  of  the  Rankine-Hugoniot  relations  in  the  viscous-layer  model 

presupposes  that  the  shock  thickness  and  the  transport  effects  behind  the  shock 

38  39 

are  not  important.  Germain  and  Ghiraud,  and  also  Chow  and  Ting,  noted, 
however,  that,  in  order  to  claim  superiority  over  the  vorticity  and  other 
second-order  boundary-layer  theories,  theories  based  on  the  viscous-layer 
model  should  also  include  the  shock  thickness  effects.  These  effects,  though 
small  for  a  high  shock  compression  ratio,  have  been  analyzed  as  corrections 
in  Refs.  38  and  39. 

A  more  recent  analysis  of  the  stagnation  region  in  the  merged-layer 

40 

regimes  is  given  by  Levinsky  and  Yoshihara.  They,  unlike  Probstein  and 
Kemp^  who  assume  a  constant  density  for  the  shock  layer  and  treat  the  shock- 
transition  zone  by  an  approximate  method,  applied  the  compressible  viscous- 
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flow  equations  along  the  axis  of  symmetry  from  the  body  upstream  to  infinity, 
including  the  shock-transition  zone.  Although  a  great  number  of  the  higher  - 
order  terms  from  the  Navier-Stokes  equations  are  retained  in  their  analysis, 
the  reduced  system  of  equations  is  self  consistent  only  to  the  leading  approx¬ 
imation.  This  is  apparent  from  the  assumption  of  a  density  field  having 
spherical  symmetry.  On  account  of  this  degree  of  approximation,  the  analysis 
cannot  be  applied  to  the  entire  domain  of  the  "fully-merged-layer  regime"  as 
originally  defined  by  Probstein  and  Kemp^  (which  permits  a  rather  thick  shock) 
without  sacrifice  of  accuracy.  The  numerical  solution  of  this  analysis  will  be 
discussed  later  in  connection  with  the  new  treatment  of  the  shock-transition  zone 
developed  in  Part  II. 

Most  analyses  described  above  have  to  make  use  in  one  way  or  another 

of  the  thin- shock-layer  assumption.  This  assumption  has  been  used  only  to 

make  the  similarity  solution  possible.  Strictly  speaking,  these  analyses  are 

no  more  accurate  than  the  leading  approximation  given  by  the  inviscid  shock- 
18-21 

layer  theory.  An  approach  making  consistent  use  of  the  thin- shock-layer 

approximation  may  therefore  describe  the  problem  with  equal  accuracy.  In 
fact,  under  this  approximation,  the  analysis  of  the  stagnation  region  embodies 
the  analyses  of  both  Sections  1  and  2. 

3.  THIN-SHOCK-LAYER  THEORY 

An  analysis  of  the  stagnation  region  based  on  a  more  consistent  application 

of  the  thin- shock-layer  approximation  has  been  given  in  a  previous  paper  by 
41 

the  author.  By  this  thin-layer  approach,  curvature,  slip  and  displacement 
all  become  negligible  higher-order  effects.  The  essential  features  of  that 
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analysis  will  be  discussed  below.  In  addition  to  reviewing  the  basic  work 
on  the  axisymmetric  problem,  previous  related  analysis  for  yawed  cylinders 
will  also  be  discussed. 


3.  1  Two  Continuum-Flow  Regimes 

In  the  thin-layer  approach,  there  are  two  important  parameters  £ 


and  (<  defined  as 

7-1 


£  = 


2-/ 


2  P  L/  a' 

K  =  e 


oo  oo  ~  fjjL  1  yU°  \ 

v  T0 


(3.1) 


The  parameter  \  ,  except  for  a  factor  weakly  dependent  on  the  reference 

temperature,  is  essentially  the  product  of  £  and  the  Reynolds  number  Reb 
used  by  Probstein  and  Kemp.  ^  For  the  analyses  of  the  stagnation  region  in 
hypersonic  flow,  the  problem  falls  into  two  regimes  (see  Fig.  1), 


Regime  I  :  0(l)  —  £  K  <  °° 


(3.  2) 


Regime  II:  0(e)  -  £  K  *  ^  0(l) 

Note  that  the  lower  limit  of  Regime  II  includes  a  major  portion  of  the  fully 
merged-layer  regime  where  the  shock  layer  and  the  shock-transition  zone 
are  comparable  in  thickness.  Also,  because  of  the  overlap  at  E  K  Z—  0(1)  , 
the  viscous-layer  regime  is  covered  by  both  Regime  I  and  II. 

3.  2  The  Higher  Reynolds  Number  Regime 

The  problem  of  the  stagnation  region  with  axisymmetry  is  reducible 
in  Regime  I,  which  includes  also  the  viscous-layer  regime,  i.e.  &K*—0(l) 


In  terms  of  £  and  K*  ,  the  domain  of  the  incipient-merged-layer  and 
fully-merged-layer  regimes  defined  originally  by  Probstein  and  Kemp  are, 
respectively,  0  (YeT )  ^  6  K  2  —  0(1 )  and  0(e^z)  (=  E  K  2  —  0(/£~) 
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This 


virtually  to  the  vorticity-inte  raction  theory  of  Hayes  and  Probstein. 
is  possible  because  under  the  thin- shock-layer  approximation,  which  is  subject 
to  an  error  of  order  &  ,  the  asymptotic  solutions  to  the  stagnation-point 

boundary  layer  also  describe  the  inviscid  (vortical)  portion  of  the  shock  layer. 
Since  the  Reynolds  number  in  this  regime  is  still  quite  high,  one  may  replace 
the  shock  boundary  by  in  the  boundary-layer  coordinate.  A  more 

detailed  discussion  of  this  reduction  is  given  in  Ref.  41  . 

The  system  of  equations  governing  the  axisymmetric  stagnation  region 
in  Regime  I  i s 

f'"  +  ff"  =  \  f  '*-  ^ 
e"+  pr?  @'  =  o 
f(0)  =  f'(O)=  9(0)  -o 
@(oa)~  I  ;  f'(co)=  Cl* 

where  the  parameter 

,  8& 

yr(3$*K  ,34) 

is  basically  the  same  as  I~L  of  Hayes  and  Probstein  (if  the  reference  temper¬ 
ature  T*  is  equal  to  the  wall  temperature  Pw  ,  can  be  reduced 

essentially  to  j(q_  ).  One  recalls  that,  in  arriving  at  Eq.  (3.  3)  or  its 
equivalent,  Hayes  and  Probstein  assume  jTL  small  and  have  to  invoke  the 
principle  that  the  flow  properties  along  the  outer  streamlines  are  not  affected 
by  the  boundary-layer  displacement.  Under  the  present  thin- shock-layer 
approximation,  this  principle  is  fully  justified,  and  the  current  controversy 
regarding  the  vorticity-displacement  effect  is  by-passed.  In  the  present 
formulation,  the  magnitude  of  ft*  is  not  necessarily  small,  in  fact,  in  the 
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viscous-layer  regime,  _Q_X  =  o(s~''4^  In  the  case  of  plane  flows,  the 
external  vorticity  vanishes  at  the  surface  and  has  relatively  small  influence 
on  the  boundary  layer.  The  treatment  of  the  viscous -layer  regime  based  on 
Eq.  (3.  3),  in  addition  to  having  removed  the  need  for  simultaneously  determining 
the  shock  stand-off  distance  with  other  unknowns,  has  the  merit  of  combining 
the  specific -heat  ratio  y  and  the  Reynolds  number  into  the  parameter  £)_  , 

3,3  The  Lower  Reynolds  Number  Regime 

Although  the  slip  and  temperature- jump  at  the  surface  may  eventually 
become  larger  at  lower  Reynolds  number,  they  cannot  exceed  the  order 

,  i.  e.  i  so  long  as  the  thin- shock-layer  approximation 


e  -j1-  CH 
•  a 


42  43 

applies.  ’  For  a  cold  body,  these  effects  may  therefore  be  neglected,  or 

separately  treated.  On  the  other  hand,  the  transport  processes  immediately 

behind  the  curved  bow  shock  give  rise  to  substantial  changes  in  the  tangential 

velocity  component  and  total  enthalpy  across  the  shock.  These  effects  resemble 

the  slip  temperature  jump  phenomena.  Probstein  in  Refs.  10  and  11  emphasizes 

the  importance  of  the  shock-thickening  effects  is,  the  merged-layer  regimes. 

41  42 

The  study  under  the  thin- shock-layer  approximation  ’  reveals,  however, 
that  the  contribution  of  the  finite  shock  thickness  is  only  secondary  and  may 
be  neglected  as  compared  to  the  more  important  "slip  effects"  at  the  shock. 

(This  will  be  further  confirmed  by  the  development  in  Section  5.3)  In  fact, 
this  is  borne  out  by  Probstein  and  Kemp's*  °  analysis  of  the  shock-transition 
zone  (although  it  was  not  clearly  brought  out  in  that  work).  Thus,  even  in 
the  merged-layer  regime  the  Rankine-Hugoniot  relations  need  only  be  modified 


The  slip-like  transport  effects  at  the  shock  have  recently  been  studied  by 
Pan  and  Probstein  in  the  leading-edge  region  of  a  sharp  flat  plate.  43 
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to  account  for  the  transport  effects  behind  the  shock,  treating  it  as  a 
surface  of  discontinuity.  Contributing  to  this  result  are  two  factors: 
first,  the  shock  thickness  in  the  regimes  considered  is  still  thin  in  com¬ 
parison  with  the  shock's  radius  of  curvature,  and  second,  the  average 
density  within  the  shock-transition  zone  is  actually  much  closer  to  the  low 
free -stream  density  than  that  in  the  thin- shock  layer  behind  it. 

For  the  axisymmetric  stagnation  region  and  a  linear  viscosity-temper ature 
relation,  an  analytical  development  is  possible  in  this  regime.  In  terms  of  the 
variable  0  =  /p  U  V  ?.Z  >  related  to  the  stream  function,  the  leading 

terms  in  the  development  for  the  tangential  velocity  component  and  the  enthalpy 
function  are  simply 


U  = 


U. 


Uoo  M>»/S 


=  K‘ 


-! 


S/2 


(3.  5) 


0  =  =r,(^y  [r,C£>i'vV 


where  F,  (t)  is  the  incomplete  gamma  function  of  order  1/3  in  the  argument 
~t  ,  and  the  parameter  Jfc,  is  equal  to  2  R^/3  ~/\  +  - /J  .  The  second 

terms  in  the  development,  which  are  proportional  to  E  ,  have  been  deter¬ 
mined  and  given  in  Ref.  41.  Except  for  a  constant  of  proportionality,  the 

velocity  ratio  (/  of  Eq.  (3.7)  is  identified  in  form  with  the  corresponding 

19,  20 


solution  of  inviscid  Newtonian  shock-layer  theory  for  £■ 


■0 


Because  of  the  dependence  on  (<"  ,  the  magnitudes  of  the  correction  terms 

are  actually  not  given  by  0(e)  but  are  considerably  larger  near  the  body 
(See  Ref.  42).  By  a  systematic  expansion  in  terms  of  y ~  *  ,  Shidlovsky 


has  recently  obtained  an  analytic  development  similar  to  that  described  above. 
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3.4  Stagnation  Region  of  a  Cylinder  at  Yaw 


The  same  thin-layer  model  has  been  applied  to  analyze  the  stagnation 
region  of  a  cylinder  at  yaw.  ^ 

In  the  case  of  an  unyawed  cylinder,  the  heat-transfer  rate  determined 
from  the  analysis  (for  £  =0.10,  Tw^/T0 — >0  )  shows  little  departure  from 

the  rate  based  on  the  boundary-layer  theory  for  «  *  as  small  as  1.0  even 
though  the  velocity  and  enthalpy  profiles  obtained  would  suggest  considerable 
difference  from  the  boundary-layer  theory.  Apparently,  the  (higher -order) 
vorticity  effect  and  the  slip-like  effects  behind  the  shock  tend  to  compensate 
each  other. 

The  analysis  for  the  yawed  cylinder  is  essentially  similar,  and  is 
particularly  simple  for  unit  Prandtl  number.  The  essential  parameter  govern¬ 
ing  the  yawed  cylinder  is 

« 2  A 


where  A  is  the  yaw  angle.  The  analysis  carried  out  for  unit  Prandtl 
number  shows  insensitivity  of  the  enthalpy  profile  to  yaw  at  both  high  and 
low  Reynolds  numbers.  This  result  suggests  a  yaw  independence  principle. 
Accordingly,  the  heat-transfer  rates  for  all  yaw  angles  can  be  correlated  as 


C, 


VS  X  MO  A 


/art;  A 

For  JH.  oC  T10  .  and  taking  into  account  the  dependence  of  the  reference 

.  .  2,  . 

temperature  on  the  yaw  angle,  the  parameter  K  sec  A  may  be  replaced 


-I 


This  rule  of  correlation  is  consistent  with  the  previous  results  obtained  by 
8  Z  8  3 

Reshotko  et  al  ’  for  the  boundary-layer  limit.  The  results  of  Ref.  42 
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■V*> 


show  that  the  correlation  is  reasonably  accurate  even  for  yaw  angles  as  large 
as  60°. 

4.  COMPARISON  OF  RESULTS  FOR  ST AGNAT ION- POINT  HEAT  TRANSFER 

The  heat-transfer  rate  is  a  quantity  of  engineering  interest,  and  can 
also  be  measured  in  experiments  conducted  in  low-density  tunnels.  Some 
controversy  appears  to  have  arisen  in  recent  years  concerning  the  pre¬ 
dicted  and  measured  values  of  this  quantity  at  the  stagnation  point  of  a  sphere. 

In  the  following,  the  heat-transfer  rates  predicted  by  the  various  theories 
will  first  be  examined. 

4.  1  Theoretical  Results 

The  theoretical  results  may  be  discussed  in  two  ranges  corresponding 
to  high  and  to  low  Reynolds  numbers. 

High  Reynolds -Number  -  It  is  advantageous  to  consider  first  the 
results  obtained  from  the  thin-shock-theory ,  which  is  not  restricted  to  the 
small  vorticity  effect.  The  result  based  on  numerical  integration  of  Eq.  (3.  3) 
for  Pr  =  0.71,  give  s  (best  fit  in  the  range  of  0  —  IT  ^  /  ,  0  ^  ^-40  ) 

=  °.5°nM- o.nnl  (4-1} 

f  8L 

where  <^.  is  the  surface  heat-transfer  rate  and  the  subscript  &L  signifies 
the  boundary-layer  limit. 

27 

With  the  exception  of  the  theoretical  results  of  Ferri  and  co-workers, 

there  is  a  reasonable  consistency  between  Eq.  (4.  1)  and  the  vorticity  effects 

predicted  by  others  based  on  the  second-order  boundary  layer  theory.  Van 

1 

Dyke  in  Ref.  24,  assumes  Pr  =0.71,  "R/T0  -  0.  20,  and  7  -  1 . 40.  In 
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27 


the  form  of  Eq.  (3.  5),  his  results  (excluding  the  temperature-jump  effect) 

2  6 

give  the  value  0.46  as  compared  to  the  coefficient  0.50.  Maslen,  and 

more  recently  Probstein  (private  communication),  consider  pr  =  1  and 

JV/]~0 — >0,  and  obtain  0.45  and  0.44,  respectively,  Ferri  and  co-workers 

give  their  results  in  terms  of  a  parameter  {b/7SARe  F  which  may  be  closely 

identified  with  .  However,  the  vorticity  effect  predicted  by  them  appears 

to  be  30  to  50%  larger  than  the  above  results. 

There  is  an  apparent  inconsistency  between  Eq.  (3.  5)  and  the  earlier 

4 

result  of  Hayes  and  Probstein  who  gave  a  coefficient  of  0.19  instead  of 
0.  50.  Their  calculation  was  made  for  Fr  =0.71,  and  a  low  value  of  /["„  • 
Instead  of  the  linear  viscosity-temper ature  relation,  as  used  in  most  other 
analyses,  the  Sutherland  law  is  used.  Clearly,  the  discrepancy  results  from 
the  difference  in  the  viscosity  law.  In  order  to  compare  the  two  results,  one 
must  use  the  appropriate  reference  temperature  in  the  linear  viscosity  law, 
which  in  this  case  can  be  taken,  according  to  the  study  of  Ref.  45  as 

(4.1a) 


' w 


T0 


*  Z  2 

and  one  must  remember  that,  implicit  in  Hayes  and  Probstein's  parameter 

_Q_  of  Eq.  (1.2),  the  wall  temperature  is  taken  as  a  reference  temperature. 

Allowing  for  a  low  value  of  Tw^~fo  ,  the  relatively  small  numerical  coefficient 

42 

given  by  Hayes  and  Probstein  can  actually  be  accounted  for.  The  apparent 
discrepancy  should  not  be  taken,  however,  to  indicate  inadequacy  of  the  linear 


In  this  example  with  a  cold  surface,  the  effects  associated  with  curvatures 
and  displacement  are  very  small.  The  ambiguity  resulting  from  the  arbitrary 
division  of  the  vorticity  and  displacement  effects  mentioned  in  Section  1  does 
not  arise . 
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representation  of  the  viscosity-temperature  relation.  Examples  in  Ref.  45 
show  in  fact  that  with  the  proper  reference  temperature  the  linear  relation 
is  adequate  for  predicting  the  heat  transfer  rate,  skin  friction,  shock  stand¬ 
off  distance  and  other  quantities  of  interest.  Rather,  the  discrepancy 
reflects  the  inappropriateness  of  using  the  wall  temperature  as  a  reference 
condition  for  the  definition  of  S\_  In  fact,  Hayes  and  Probstein's  results 


lead  to  an  infinite  rate  of  heat  transfer  in  the  limit  of 


TL 


T0 


-+0  • 


Low  Reynolds  Number  -  Heat-transfer  results  pertaining  to  the  merged- 

layer  regimes  have  been  obtained  by  Probstein  and  Kemp,  ^  Levinsky  and 
41  42 

Yoshihara,  and  the  author.  For  these  regimes,  the  thin- shock-layer 


CH  =  3/fc 


(4.2) 


approach  yields  the  first  approximation  consistent  with  Eq.  (3.  5)  as 

E] 

which  approaches  the  free -molecule  limit  corresponding  to  a  unit  accommoda¬ 
tion  coefficient.  One  notes  that  the  effect  of  the  specific  heat  ratio  7  will 
only  appear  as  a  correction  to  the  value  of  oH  given  above,  which  is  a  function 
only  of  ■  Therefore,  flows  with  different  values  of  7  should  correlate 
reasonably  well  in  terms  of  the  parameter  |<2  in  the  merged-layer  as  well  as 
in  the  viscous-layer  regimes  (see  the  numerical  solutions  given  in  Ref.  41  and 
42).  In  fact,  the  heat- transfer  rates  based  both  on  the  solutions  of  Levinsky 


and  Yoshihara  (  y  ;  -— )  and  that  of  Probstein  and  Kemp  (constant-density 

model  with  7 -  ),  agree  reasonably  well  with  Eq.  (4.  2).  The  corresponding 

shock  stand-off  distance  (to  be  more  precise,  the  thickness  of  the  shock  layer) 

predicted  by  Probstein  and  Kemp  is,  howevei  ,  much  higher  than  that  obtained 

41 

by  the  present  approach,  primarily  because  of  their  constant-density 
assumption 
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4,2  Experimental  Re  suits 

Measurements  of  stagnation-point  heat  transfer  have  been  made  on 

46 

spheres  in  hypersonic  rarefied  gas  flow  by  Neice  and  co-worker,  Ferri 

27  47  48  49  50  51 

and  co-workers,  ’  ’  Wittliff  and  Wilson,  ’  Varig,  Bloxsom  and 

52 

Rhodes,  and  Hickman  and  Giedt  (unpublished  but  reported  in  Ref.  12). 

Similar  data  for  yawed  and  unyawed  cylinders  were  also  obtained  by  Vidal 
50 

and  Wittliff.  While  discrepancies  do  exist  among  the  data  reported,  as 
well  as  among  analyses,  the  general  trend  of  the  data  does  bear  out  the 
theoretical  predictions.  The  disagreement  between  theory  and  experiment 
is  found  mainly  in  the  vor ticity-interaction  regime  where  the  departure 
from  the  boundary-layer  prediction  is  actually  quite  small. 

27  47 

In  Fig.  4.  1,  the  data  for  spheres  of  Ferri  and  co-workers  ’ 
which  encompasses  both  the  high  and  low  Reynolds  number  regimes,  are 
reproduced  along  with  theoretical  predictions  based  on  their  analysis  of 
Ref.  27.  The  data  were  taken  in  tests  with  a  stagnation  temperature 
T0  =  1  280°K  and  test  flow  Mach  numbers  of  5  and  8.  If  one  assumes 


vibrational  equilibrium,  the  value  of  t  behind  the  shock  should  be  about 


0.13  to  0.14.  The  expe  rimental  data  are  pre  sente  d  as  p/p  vs  Re  F  , 
where  is  the  stagnation- point  heat-transfer  rate  measured  on  a  small 
sphere  and  p  is  the  heat-transfer  rate  measured  on  a  considerably 
larger  sphere  corresponding  to  vanishing  vorticity.  The  Reynolds  number 


ReF  may  be  related  to  |<  as 
Op  =  , 


_  )  (Ta  M*\  K 

fz  £  /TV  Th  ‘  Mo  6* 


(4.  3) 
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Included  in  the  same  figure  for  comparison  is  the  theoretical  curve  based 

on  the  thin- shock-layer  approach  described  in  the  last  section  (for  both 

Regimes  I  and  II)  calculated  for  Pr  =0.71,  0,  and  &  =0,13 

Also  included  in  the  figure  are  the  relatively  recent  measurements  of  Vidal 
50 

and  Wittliff,  Their  experiments,  performed  in  a  hypersonic  shock  tunnel, 
extend  the  measurement  ot  a  range  corresponding  to  the  fully-mergcd-layer 
regime  and  beyond.  Although  the  shock-tunnel  data  of  Vidal  and  Wittliff  are 
obtained  for  a  test  section  stagnation  temperature  as  high  as  T0  =  2500  — 
3000°K,  chemical  reaction  is  expected  to  be  frozen  both  in  the  shock-layer 
as  well  as  on,.the  (noncatalytic)  surface.  There  was  nevertheless  a  consid¬ 
erable  concentration  of  atomic  species  present  in  the  test  flow  of  the  experi¬ 
ments  due  to  freezing  in  the  nozzle  expansion.  In  the  data  correlation  of 
Ref.  50  the  energy  associated  with  the  frozen  atomic  species  has  been 
subtracted  from  the  total  enthalpy. 

The  general  trend  of  agreement  of  these  data  with  theory  is  quite 
evident  from  Fig .  4.1.  In  the  higher  Reynolds  number  range  (  Re  —  500), 
there  is  a  rise  over  the  boundary-layer  value  due  to  the  external-vorticity 
effect.  In  the  lower  Reynolds  number  range  (  500),  the  heat  transfer 

ratio  <2.  /j .  begins  to  decrease  and  soon  goes  below  unity,  approaching 


u/  -eL. 

the  free  molecule  limit.  If  one  considers  the  scatter  of  the  data  and  the 
degree  of  approximation  in  the  analysis,  the  agreement  is  reasonable.  To¬ 
ward  the  high  Reynolds  number  end,  the  data  of  Ferri  and  co-workers  agree 


Only  data  in  the  range  of0.10£f<'1f=l  of  Ref.  50  are  included  in  Fig. 
4.  1  . 
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well  with  their  own  theoretical  prediction,  but  appear  to  be  higher  than  that 
predicted  by  the  thin- shock-layer  approach  (as  well  as  the  others).  As  far 
as  the  fractional  increase  in  heat  transfer  (t-f «.)/*«.  is  concerned, 

measurements  of  Ferri  and  co-workers  are  consistently  higher  than  the 
writer's  prediction  by  30  to  50%,  The  unsolved  discrepancy  suggests  that 
effects  of  vibrational  (and  perhaps  also  rotational)  relaxation  in  the  experiments 
on  the  vorticity  effect  should  be  examined. 

Heat-transfer  data  for  sphere  obtained  by  Hickman  and  Giedt  (reported 
by  Van  Dyke)  at  Mach  numbers  from  2  to  6  in  a  low-density  wind  tunnel  have 
been  compared  with  the  second-order  boundary  layer  analysis  by  Van  Dyke.  ^ 

The  measured  departure  from  the  boundary-layer  theory  is  seen  to  exceed  the 
prediction  by  a  factor  of  two  or  more.  It  has  been  noted  in  Section  1  that 
the  slip  and  temperature  jump  effects  in  Van  Dyke's  calculation  are  comparatively 
large.  As  indicated  by  Van  Dyke  himself,  simply  deleting  them  would  bring 
better  agreement  with  experiment. 

Heat-transfer  rates  along  the  stagnation  streamlines  of  yawed  and 

unyawed  cylinders  have  bean  measured  in  the  low-density  hypersonic  shock 

45 

tunnel  and  compared  with  the  analysis  of  Cheng  and  Chang  by  Wittliff 
50 

and  Vidal.  The  comparison  bears  out  the  validity  of  the  independence 
principle  mentioned  in  Section  3.  3, 
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PART  II.  EXTENSION  OF  THE  THIN  SHOCK  LAYER  APPROACH 


The  study  in  Part  I  reveals  clearly  that  the  thin- shock-layer  approxima¬ 
tion,  appropriately  modified  to  take  account  of  the  transport  effects,  simplifies 
and  unifies  analyses  in  both  the  low  and  high  Reynolds  number  regimes.  In 
Part  II,  this  approach  will  be  extended  to  study  flow  regions  other  than,  or  re¬ 
moved  from,  the  stagnation  region;  a  treatment  of  the  shock-transition  zone  is 
included.  In  Section  5  below,  the  basic  flow  model  will  be  developed,  the  prob¬ 
lem  formulated,  and  the  adequacy  as  well  as  limitations  of  the  approach  dis¬ 
cussed.  Section  6  will  present  and  discuss  a  numerical  method  based  on  a  finite 
difference  approximation  for  solving  the  shock-layer  equations.  The  method  is 
applied  in  Section  7  to  study  flows  over  wedge,  paraboloid  and  other  simple 
nonslender  shapes.  Finally,  in  Section  8,  the  thin- shock-layer  approach  is 
generalized  to  include  nonequilibrium  dissociation  and  specific  analysis  is  made 
of  the  stagnation  region. 

5.  THE  THIN  SHOCK  LAYER  APPROACH  TO  THE  VISCOUS  BLUNT  BODY 
PROBLEM 

In  this  section,  the  basic  assumptions  and  the  flow  models  underlying  the 
thin-layer  approach  will  be  discussed;  equations  governing  the  shock  layer  and 
the  shock-transition  zone  will  be  derived;  and  the  initial-value  problem  formu¬ 
lated.  As  examples  to  demonstrate  adequacy  of  the  general  formulation,  solu- 

12  53 

tions  are  tested  against  an  exact  numerical  solution  of  Van  Dyke  and  Swigart 

for  a  paraboloidal  shock  in  the  inviscid  limit,  and  against  the  numerical  solution 

40 

of  the  stagnation  region  obtained  by  Levinsky  and  Yoshihara  in  the  merged- 
layer  regime. 
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5.  1  The  Two-Layer  Flow  Model 


For  the  purpose  of  the  present  analysis,  the  flow  field  around  the  body 
may  be  divided  into  an  inner  and  an  outer  region  as  illustrated  in  Fig.  5.  1: 
the  shock  layer  and  the  shock-transition  zone.  For  convenience,  the  surface 
separating  the  two  regions  will  be  referred  to  as  the  "shock  interface".  The 
compression  ratio  o/p  across  the  shock-tr ansition  zone  is  assumed  to  be 
high.  Both  regions  will  be  assumed  to  be  thin  in  comparison  with  a  typical 
dimension  of  the  body  (say  the  nose  radius  or  a  lateral  dimension);  this  assump¬ 
tion  is  consistent  with  high  compression  ratio  assumed. 

With  the  exception  of  the  study  of  Section  7  which  deals  with  the  nonequi¬ 
librium  dissociation,  an  ideal  gas  with  constant  specific  heats  is  assumed. 

In  the  subsequent  analyses,  a  linear  viscosity-temperature  relation  may  be 
adopted.  The  formulation  presented  applies  to  plane  and  axisymmetric  steady 
flows  around  smooth,  nonslender  bodies.  Applications  are  made  mainly  to 
highly  cooled  surfaces  with  low  wall-to- stagnation  temperature  ratios. 


Let  4,  denote  the  distance  along  the  body  surface  and  y-  the  distanc* 
along  the  outward  normal  from  the  surface;  <£  and  form  a  pair  of 

orthogonal  curvilinear  coordinates  (refer  Fig.  5.  1).  It  is  assumed  that  both 
the  shock  layer  and  the  shock- transition  region  are  thin  so  that  -g— -  2.C  -g— . 


Implicit  in  the  assumption  of  constant  specific  heats  is  the  assumption  that 
internal  degrees  of  excitation,  such  as  molecular  vibration  and  dissociation 
do  not  undergo  the  processes  of  relaxation.  Under  this  condition,  one  may 
assume  a  zero  '  bulk  viscosity". 

The  nonslender  bodies  admissible  for  the  present  study  can  be  pointed,  but 
the  surface  angle  must  not  be  small  and  (except  the  apex  of  a  pointed  body)  the 
surface  curvature  must  be  finite. 
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The  Navier-Stoke s  equations  governing  the  compressible,  viscous,  heat- 
4  24 

conducting  flow  ’  may  then  be  reduced  to  the  following  form,  neglecting 
the  transverse  and  longitudinal  curvature  effects  and  retaining  only  the 
highest-order  derivatives  with  respect  to  <y.  in  the  transport  terms. 

=  o  ^ 

^t6cu-z  +■  e  f 

^  +  ^  2  /'fU  2  +  3  VyJ  Y  / 

Here  U.  is  the  velocity  component  along  ^  \  /\r  ,  velocity  component 
along  ^  ;  p  ,  the  density,  -p  ,  the  pressure;  X  ,  the  specific  enthalpy; 

/<  ,  the  coefficient  of  viscosity;  2  the  distance  of  the  body  surface  from 

/  fi 

the  axis  of  symmetry;  1 6,  ,  the  longitudinal  curvature  of  the  surface—^ —  . 
The  index  -J  is  zero  for  plane  flow  and  unity  for  axisymmetric  flow.  Pr 
is  the  Prandtl  number.  These  flow  variables  are  further  related  through  the 
equation  of  state  -p>=ppT,  A-cpA’  and  the  relation  between  yU.  and  T  . 

In  the  above  equations,  both  the  transverse  curvature  effect,  as  well  as  the 
higher -order  longitudinal  curvature  effects,  have  been  neglected. 

The  system  of  Eq.  (5.  1  )  which  has  been  given  previously  by  Hayes  and 

4 

Probstein  is  valid  in  the  shock  layer  and  in  the  shock-transition  zone  to 
within  errors  of  the  order  A  !ol,  and  S^/z  ,  respectively.  To  the  same  degree 
of  accuracy,  the  coordinates  (  <^  ,  y>  ),  and  the  corresponding  velocity  com¬ 
ponents  (  U  ,  W  ),  in  the  system  of  Eq.  (5.  1)  may  also  be  interpreted  as  the 
coordinates  and  velocities  parallel  and  normal  to  the  shock  interface  instead 


of  to  the  body  surface .  The  terms  —  p  IGC  U.  nt  in  the  second  of  Eq.  (5.1)  and 
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dnr  . 
dtf. 


the  term  U  ■— -  in  the  third  of  Eq.  (5.1)  actually  belong  to  the  higher  orders 


and  can  be  neglected.  One  may  note  that  the  system  of  Eq.  (5.1)  is  parabolic, 
To  be  more  precise,  Eq.  (5.  1)  has  a  degenerated  characteristic  pair 

—  t.  oo  ,  similar  to  the  boundary-layer  equations  of  compressible  flow. 
Equations  (5.1)  can,  however,  be  simplified  further  in  the  shock-layer  region 
and  in  the  shock- transition  zone. 

5.  2  Governing  Equations  in  the  Shock  Layer 

Because  of  the  high  compression  ratio  across  the  shock-transition  zone, 
one  has  in  the  shock  layer,  in  addition  to  the  simplification  — 2.2, 


Poo/e 


^  I 


7 


(5.2) 


This  implies  that  the  normal  velocity  component  within  the  shock  layer  is  small, 
or  more  specifically, 


'ir 

Uco 


(5.2a) 


One  notes  that  the  ratio  /Vy  U  is  usually  small  in  the  shock  layer,  but  there  is 

no  requirement  that  it  has  to  remain  small.  This  observation  is  essential  with 

respect  to  the  validity  of  the  following  equations  in  the  stagnation  region.  Anti- 

?.  .  z 

cipating  that  P'  is  of  the  order  of  jO  U^Abtb  j$  ,  and  U  is  of  the  order 

and  less  within  the  shock  layer,  Eq.  (5.  1)  can  be  further  simplified  to 


If  terms  of  order  ^A+ds)/d  in  the  Navier-Stokes  equations  were  also  retained, 
the  partial  differential  equations  would  contain  such  terms  as  U.Xy  .  'V'Ky  •  etc  • 
and  would  therefore  appear  to  be  hyperbolic.  One  still  has  an  initial- value  prob¬ 
lem,  however. 

The  subsequent  formulations  of  the  problem  in  the  two  flow  regions  are  com¬ 
parable  to  solving  Eq.  (5.  1 )  by  a  procedure  of  inner  and  outer  expansion  under 
the  assumption  of  high  compression  ratio  across  the  shock. 
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(P  u  *  (p*  l\  =  0 

r*+  ^)u  = 

fy  *  *cP“-Z  =° 

p(u  s|)H  =jfiT  [n-<i-Pr)-§-^ 

where  H  is  the  total  enthalpy.  The  system  Eq.  (5.  3)  which  is  to  be  applied 

in  the  shock  layer  is  subject  to  a  fractional  error  of  the  order  A/a-  ,  or  £  . 

This  system  differs  from  the  boundary-layer  equations  only  in  one  respect, 

namely,  the  pressure  variation  across  the  layer  is  no  longer  negligible  and  is, 

in  fact,  accounted  for  here  by  the  same  equation  as  for  the  inviscid  shock  layer. 

One  may  note  that  in  the  present  formulation  the  tangential  pressure  gradient 

aty.  ,  which  has  been  generally  regarded  as  a  higher-order  term  in  the  exist- 

41 

ing  inviscid  theory,  is  retained.  As  has  been  observed,  this  term  will  be 
important  near  the  body  surface  when  the  Reynolds  number  is  high,  since  the 
flow  speed  near  the  base  of  the  shock  layer  is  of  the  order  It  instead  of 
according  to  inviscid  theory.  The  term  -fy  is  therefore  retained  as  a 
principal  term  in  order  to  insure  uniform  validity  of  the  equation  at  the  higher 
Reynolds  numbers.  This  is  essential  for  a  correct  description  of  the  boundary- 
layer  phenomena  under  the  present  formulation. 

Before  turning  to  the  shock-transition  zone,  one  should  recall  that  there 
are  two  well  recognized  shortcomings  of  the  inviscid  shock-layer  theory  which 
would  also  appear  to  affect  the  present  approach.  One  is  the  breakdown  of  the 
approximation  at  a  certain  critical  region  on  a  convex  body,  corresponding  to 
the  shock-layer  "separation"  at  the  zero-pressure  point.  ^  This  difficulty 


AF-1285-A-10 


29 


is  not  found,  however,  with  certain  shapes  such  as  the  cones,  paraboloids  and 

hyperboloids  to  be  studied.  In  any  case,  the  viscous  effect  tends  to  reduce  the 

flow  speed,  thus  the  centrifugal  action,  and  delay  the  separation.  The  second 

difficulty  has  to  do  with  the  slow  convergence  of  the  Newtonian  series,  familiar 

54 

from  the  work  of  Chester,  which  renders  the  leading  approximation  worthless. 

This  shortcoming,  however,  is  not  applicable  to  the  present  formulation.  As 

4  55 

is  apparent  from  a  discussion  by  Van  Dyke,  ’  the  inaccuracy  may  be  traced 
back  to  the  tangential-pressure-gradient  term  disregarded  in  the  leading 

approximation.  The  present  formulation,  which  has  included  as  a  principal 
term,  averts  this  difficulty  as  will  be  demonstrated  by  subsequent  examples  in 
Section  5.6. 


5.  3  Governing  Equations  in  the  Shock-Transition  Zone 

From  the  existing  analyses  of  the  shock-wave  structure  in  one  dimen- 
55-57 

sion,  one  may  anticipate  the  normal  velocity  component  to  vary  more 


or  less  in  an  anti- symmetric  manner  with  respect  to  y- 


Detween  tne  two  limits 


0(t>  -  u^/s  -  1  (5'4) 

It  follows  that,  excluding  a  layer  of  thickness  £  5  near  the  shock  interface 
(see  Fig.  5.1),  one  has  in  the  shock-transition  zone 

-  0(1)  (5.5) 


Anticipating  pnr  ~  p  U^^yrvp  ,  it  also  follows  from  Eq.  (5.4)  that  in  the 
shock-transition  zone  (again  excluding  the  region  near  the  shock  interface) 


00 


(5.6) 


That  is,  the  density  in  the  shock-transition  zone  is  generally  at  the  (low)  free- 
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stream  level.  One  must  note  the  fact  that  an  anti- symmetric  normal-velocity 
profile  will  give  a  highly  asymmetric  density  profile  for  high  compression  ratio, 
since  p  cQ.  ~ •  This  explains  why  the  average  density  through  the  shock- 
transition  zone  is  so  low  compared  to  the  density  at  the  interface.  In  any  case, 
the  assumptions  related  to  Eqs.  (5.  5)  and  (5.  6)  are  readily  verified  a  posteriori. 

For  the  study  of  the  shock-transition  zone,  it  will  be  more  convenient  to 
employ  the  system  of  orthogonal  coordinates  in  which  the  surface  ^  =  0 
coincides  with  the  shock  interface. 


Using  Eqs.  (5.  5)  and  (5  6)  and  the  thin-layer  assumption,  the  system  of 
Eq.  (5.1)  is  then  reduced  simply  to 


(P^)f  =  0 

PAruf  =(>  Uy) 


+  P'^y  =  T  O 


(5.7) 


which  can  be  simply  integrated  to  yield  a  system  of  first-order  ordinary  differ¬ 
ential  equations  governing  the  shock-transition  zone 

P^  '  Pi**!  ^ 


■p'  +  p,  AT,  — f  M  'Vy  -  P,  or  \ 

«  (*  *  --z—)  i-  *  ^  J 

where  the  subscript  1  refers  to  the  upstream  condition,  with  p  s  p  ,  A,  =  ^ 
U.,  =  and  Arf  =  —U^ydi^S. 

Noting  that,  immediately  behind  the  shock,  the  system  of 

equations  Eq.  (5.8)  after  neglecting  terms  of  order  of'  is  seen  to  include  the 


modified  Rankine -Hugoniot  shock  relations  of  Refs.  41  and  42: 
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Pl^Z  =  P,^, 

^2  Z  P,^ 

P,%  =  0“  “f)2  ^  (5'9) 

P/^  (Hz  -H,)  -\j5~  H  -('-&)  t]  | 

where  the  subscript  2  refers  to  the  condition  at  the  shock  interface. 

If  the  alternative  coordinate  system  is  used  in  which  the  surface  y-  -  0 
coincides  with  the  body  instead  of  with  the  interface,  the  first  of  the  system  of 


Eq.  (5. 7)  should  be  written  as 


=  0.  The  term 


dA 

dtf. 


u 


accounts  for  the  fact  that  maximum  gradients  do  not  occur  in  the  same  direction 
as  the  normal  to  the  body- surface,  and  is  important  near  the  interface  where 
or  becomes  small  and  of  order  £  Ugg  .  The  appearance  of  this  term  does  not 
materially  complicate  the  matter.  The  whole  system  of  the  governing  equa¬ 
tions  can  be  written,  without  increasing  the  degree  of  inaccuracy,  in  a  form 
identical  to  Eq.  (5.  7)  with  sir  replaced  by  ‘ 


/  _  dA  \ 
a*  7 


The  simplified  Eq.  (5.  7)  is  not  strictly  valid  in  the  neighborhood  of  the 
shock  interface,  i.  e.  ,  in  the  shaded  region  of  Fig.  5.1.  However,  the  equa¬ 
tions  after  integration  with  respect  to  y-  (the  distance  along  the  normal  to 
the  interface),  that  is  Eqs.  (5.8)  as  well  as  shock  relations  Eqs.  (5.9),  are 
valid  uniformly.  The  remainder  in  each  of  these  integrated  equations  causes 
an  error  in  the  solution,  which  is  at  most  of  the  order 


■7 


A. 

CL 


ju  -4- 


(5.  10) 


'  Poo 

where  the  integral  is  taken  over  the  shock-transition  zone. 

One  should  note  that  the  modified  Rankine -Hugoniot  shock  relation 
Eq.  (5.9)  does  not  contain  any  information  on  the  shock-wave  structure.  A 
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similar  set  of  shock  relations  has  been  proposed  previously  by  Sedov, 

Michailova  and  Chernyi"^  for  the  case  in  which  the  transport  effects  are 
weak.  The  theoretical  basis  for  neglecting  the  shock  thickne  ss -curvature 
effects  has  not  been  provided,  however.  It  is  interesting  to  note  that,  in 
view  of  the  appearance  of  the  velocity  and  enthalpy  gradients  in  these  modi¬ 
fied  shock  relations,  the  tangential  velocity  and  the  total  enthalpy  will  change 
across  the  shock  in  a  manner  quite  similar  to  the  slip  and  temperature  jump 
across  the  Knudsen  layer  next  to  a  body  surface.  With  Eq.  (5. 9)  as  an  outer 
boundary  condition,  the  shock-layer  region  can  be  analyzed  independently  of 
the  shock-transition  zone.  Before  formulating  the  shock-layer  problem, 
however,  thetreatment  of  the  shock-transition  zone  will  be  further  simplified. 

5.4  Treatment  of  the  Shock-Transition  Zone 

The  equations  governing  the  shock-transition  zone  have  been  reduced, 
in  a  manner  consistent  with  the  thin- shock-layer  approximation,  to  a  system  of 
first-order  ordinary  differential  equations,  namely,  Eq.  (5.8).  One  notes  that 
that  system  has  satisfied  the  flow  condition  upstream  of  the  shock,  and  that  the 
boundary  conditions  for  this  first-order  system  are  to  be  furnished  by  the  shock- 
layer  solution  at  the  interface.  In  the  following,  the  shock-transition  zone 
problem  will  be  further  reduced  under  the  assumption  of  a  Prandtl  number  of 
3/4  to  the  solution  of  a  single  first-order  differential  equation. 

It  is  convenient  to  introduce  the  nondimensional  quantities 


”  ) 

nr  = 

nr 

<  ’ 

1 

i 

l 

j.  (5.11) 

-p'  - 

r 

T 

r 

\ 

| 

/  ss 

To* 

P  ”  P,  J 
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To  eliminate  the  viscosity  /A  ,  one  introduces  the  transformation 


<£, 


*)  =  T  uoo  j 

0  s  -f-  /?n  [  (  d(2-\ 


(5.12) 


J'fx 

where  gives  the  location  of  the  shock  interface.  In  terms  of  this  new  vari- 

ble  q  and  the  dimensionless  variables  U.  ,  qp  ,  etc.,  the  system  Eq.  (5.9) 
is  reduced  to 

P  =  I 


+  nr  +  nr-  - 1 


u  *  T  M 


i 

>  (5.13) 


H  +  ^pr  (h  +  (Pr  ~  ■ )  <C^A'  ^  °r  ~ 'A'-/n/  nr  "]  _  =1 


J 


where 


—  —  _  2  ,  J  ,,  _ 

H  =  T  +  /C^  jj  U  <-  y*sr»  p  sir 


To  complete  this  system  of  differential  equations,  one  requires  the  equation  of 
state  in  terms  of  the  new  variables 

PT 


*fy  "  £ 


yuy^y  * £ 


(5.14) 


Note  that  the  equations  governing  the  tangential  velocity  U  can  be  inte¬ 


grated  independently  of  the  rest  of  Eq.  (5.  13)  to  give 


U. 


=  '-('-“O' 


A  „ 

l'7i 


(5.  15) 


where  is  the  boundary  value  for  U.  to  be  furnished  by  the  shock-layer 

solution  at  the  shock  interface.  For  a  Prandtl  number  of  3/4,  having  deter¬ 
mined  U.  ,  the  last  of  Eq.  (5.  13)  governing  the  total  enthalpy  H  can  also  be 
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simply  integrated,  giving 


H  =  /-(/-  Hx)  e_5?  +  Z  a*'Z/3  (7-  £2/)  e*  (l  -  e^y 

--|-^z/5(/-Z72)'e_''1  (l-e~SnA) 


(5.16) 


where  is  the  boundary  value  for  H 

In  passing,  one  may  note  that  the  requirement  of  continuity  of  the  tangential 
velocity  U,  and  total  enthalpy  H  across  the  interface,  as  is  implicit  in  the 
provision  of  and  H ^  by  the  shock-layer  solution,  also  guarantees  the  con¬ 
tinuity  of  the  gradients  U  and  H  in  the  present  formulation. 

7  f 

To  eliminate  ^  in  the  second  of  Eq.  (5.  13),  one  makes  use  of  the  equation 
of  state  (5.  14).  With  the  help  of  the  definition  of  J-f  and  the  continuity  relation 


p  =  — -  ,  the  analysis  of  the  problem  is  finally  reduced  to 
'  at 


(l-t)0'+(c*>2f)(H-xu»'a)6u  =  I 


(5.  17) 


where  H  and  Li  have  been  determined  by  Eqs .  (5.  1  5)  and  (5.  1  6).  Note  that 
the  magnitude  of  a}'  varies  from  unity  at  upstream  infinity^ — ^-oo)  to  order 

£,  near  the  shock  interface  ^  Kj - >  dj  .  Concerning  the  term  associated  with 

in  Eq.  (5.  17),  one  must  notice  two  essential  facts.  First,  as  the  interface 
is  approached,  this  term  (being  inversely  proportional  to  at  )  becomes  the 
most  important  one  in  the  equation  and  cannot  be  dropped  even  though  it  is  pro¬ 
portional  to  £  .  Second,  with  this  term  retained,  the  condition  of  continuity 

of  the  normal  velocity  at  the  interface  is  automatically  satisfied  (to  the  leading 
approximation)  by  the  solution  to  Eq.  (5.  17).  This  follows  from  the  observation 
that,  when  at"  becomes  of  order  £  ,  Eq.  (5.17)  gives 

AT  &  t  esc? (hz  -  u* J  +  0(£2) 
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which  is  essentially  the  same  value  provided  by  the  shock- layer  solution  at 
the  shock  layer's  outer  edge. 

Taking  advantage  of  the  first-order  continuity  of  or  ,  one  may  prescribe 
^ ~fi_  at  t)  =  0  as  the  boundary  condition  for  Eq.  (5.17).  This  will  thus  insure 
continuity  (matching)  in  both  or  and  OT  within  the  degree  of  approximation 

ff 

considered. 


The  boundary  condition  for  Eq.  (5.  17)  may  now  be  written  as 


d'or  _  4  cu _  ( /*  \  {3or\  /  - 

drj  \'MoJz  \  ^  5-11 

where  and  (?or/dyjz  are  to  be  provided  by  the  solution  to  the  shock- 

layer  equations. 

It  may  be  noted,  however,  that  if  one  were  to  prescribe  or  instead  of 
cLnrfciTj at  the  shock  interface,  the  continuity  in  both  fir  and  its  first  derivative 

could  not  be  generally  achieved.  This  is  because,  on  account  of  the  term 

2  _  _ 

in  Eq.  (5.17),  even  a  change  of  order  6  in  At  at  I'j  =  0  will  critically 
affect  the  leading  approximation  of  d,fir/drj  at  f]  =  0. 

After  determining  fir  from  Eq.  (5.  17)  and  (5.  18),  p  can  be  obtained 
from  \Jfifi  >  7”  from  the  definition  of  fj  ,  and  from  the  equation  of  the 
state.  The  coordinate  O^,  corresponding  to  fj  may  finally  be  determined 


In  the  stagnation  region,  (d/lfi/dr, on  the  right-hand  side  of  Eq.  (5.  18)  can  be 
evaluated  from  (refer  to  List  of  Symbols) 

the  normal  velocity  Or  determined  from  the  shock-layer  solution  is  based  on 
a  coordinate  system  with  'U-  =  0  coincide  with  the  body  surface,  (  in 

Eq.  (5.18)  should  be  replaced  by  .  f  2 
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from  the  inverse  of  the  transform  Eq,  (5.17). 

Finally,  one  notes  that,  in  the  limit  of  infinite  Reynolds  number  the 

Rankine -Hugoniot  relations  hold,  and  — >  H z — >1  •  In  this  limit,  the 

solutions  given  by  Eqs.  (5. 15)  and  (5.16)  degenerate  to  U - I  and  H  — >-/  ; 

and  Eqs.  (5.17)  and  (5.  18)  reduce  to 

-  £  )  AT  +  ~  -  I 
'  fir 

If  =  °  (i'0) 

The  solution  for  /ir  may  then  be  obtained  by  a  separation  of  variables  and 

57  58  59 

corresponds  to  that  of  Becker  and  others.  ’  ’  In  the  general  case,  a 

simple  quadrature  of  Eqs.  (5.  17)  does  not  seem  to  be  obvious  in  view  of  the 
dependence  of  H  and  U.  on  f)  The  numerical  integration  of  Eq.  (5.17) 

is,  however,  rather  straightforward.  An  example  is  given  and  discussed  in 
Section  5.6. 

In  passing,  it  may  be  pointed  out  that  the  present  analysis  of  the  shock- 
transition  zone  differs  from  the  usual  one- dimensional  analysis  only  in  the 
boundary  condition  at  the  inner  edge.  It  provides  a  generalization  to  include 
downstream  nonuniformity  under  the  assumption  of  Sa/a.  4.4.  ( 

5.  5  Formulation  of  the  Shock- Layer  Problem 

To  complete  the  description  of  the  shock-layer  region,  one  has  to  specify 
the  boundary  conditions  at  the  body  surface.  The  velocity  slip  and  temperature 


A  "uniformly  valid  approximation"  to  the  solutiori  of  Fq.  (5.17)  is  neverthe¬ 
less  possible  in  which  Eq.  (5.17)  is  solved  with  H  and  TL  taken  as 
constants  evaluated  at  f)  =0. 
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jump  at  the  surface,  according  to  the  theories  of  Maxwell  and  Smoluchowski 
(see,  for  example,  Ref.  60),  are  proportional  to  the  local  mean-free  path  and 
the  gradients  of  velocity  and  temperature  at  the  surface.  These  effects  have 
been  investigated  by  many  within  the  framework  of  the  boundary-layer  theory 
(see  Section  1  of  Part  I).  Making  use  of  expressions  for  the  momentum  and 
energy  fluxes  provided  by  the  Navier-Stoke  s  equations,  one  arrives  at  the 
following  estimates,  in  terms  of  the  skin-friction  and  heat-transfer  coefficients 


Cf  and  CH 

T-Tv 


T. 


w. 


Cp, 


t°w  7rRTyj  zT 

,  3«  j 


/*  *3 


(5.20) 


i£L 


ye 

'  I  o 


U; 


•  C 


^ o  '  ^  r° 

The  magnitudes  of  C H  and  c„  are  small  in  most  flow  regimes  and  cannot  ex¬ 
ceed  unit  order.  Hence,  for  a  thin  shock  layer  {  £  46  !  )  and  a  cold  surface 

,  the  velocity  slip  and  temperature  jump  at  the  surface  are  small. 
In  the  present  study,  non- slip  boundary  conditions  will  be  used.  A  correction 
for  these  effects  can  be  made  afterwards,  if  necessary. 

For  the  present  approach,  it  is  advantageous  to  use  the  von  Mises  coor- 


dinates  (  ),  where  V  is  the  stream  function  defined  by 


3£ 


=  -(Z*r  Z)  p*r, 


■IT.  ,  (W)V 
dr 


(5.21) 


The  equations  governing  the  thin  shock  layer,  Eq.  (5.  3),  may  then  be  recast 
into  the  form 


CrU 


djv  _ _ 

3Y  ~  (%irl) 


uiJT -h'rrZ?'’  «  -2-(MpVL±'LS-  ±  2±L 
u  d*  ar  \/*r  3T/  p  s*. 


(5.22) 


9H  _  /„  7 T')  3_ 
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The  inner  boundary  condition  is,  at  -  0, 


u  -  0  , 


H  -  H 


w 


(5.23) 


and  the  outer  boundary  condition  is  provided  by  the  modified  Rankine -Hugoniot 
shock  relations,  Eq.  (5.  9),  as 
P  =  Poo  5 


(5.24) 


“  =  -(Zvlj 

These  equations  are  to  be  applied  at  the  shock  interface  where  IT—  ^  ifrif. 
Note  that  under  the  present  approximation,  2  can  be  taken  as  either  the  shock 
surface  or  the  body  surface.  The  system  of  Eqs.  (5.22),  (5.23)  and  (5.24)  com¬ 
pletely  determines  the  flow  quantities  U  ,  H  and  -p/  in  terms  of  the  variables 
<)£  and  for  a  given  surface  Z.  -  Z  ('%)  .  The  normal  velocity  /ir  ,  the 
coordinate  ■y-  corresponding  to  V  ,  as  well  as  the  location  of  the  shock  inter¬ 
face,  can  be  determined  through  Eq.  (5.21). 

Of  interest  is  the  behavior  of  the  solution  to  Eqs.  (5.  22)  to  (5.  24)  in  the 
formal  limit  of  zero  Reynolds  number.  On  account  of  the  predominance  of  the 
transport  effects  at  this  limit,  solutions  to  Eq.  (5.22)  will  be  characterized  by 

a  uniform  distribution  of  /y'  ,  T  ,  -p*  ,  ^  <  and  K>  it~  across  the  shock 

o  * 

layer.  Hence,  the  values  of  the  last  three  quantities  at  the  body  surface  will  be 
provided  by  the  shock  relations  Eq.  (5.24),  giving 


■Pw  =  1**  =  Px  £  ^  > 
v  =  CM  |~)2  =  Pco  ul 
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These  limiting  values  are  no  more  nor  less  than  the  pressure,  stress  and  heat- 
transfer  rate  in  a  hypersonic,  free-molecule  flow  for  unit  surface  accommoda¬ 
tion  coefficients. 

The  errors  introduced  in  the  application  of  Eqs.  (5.22)  to  (5.24)  belong, 
respectively,  to  the  orders  A/t 7.  ,  and  8Z>V  t/a.  .  In  the 

present  study,  only  cold  surfaces  will  be  considered.  The  analyses  will  be 
further  limited  to  regimes  where  the  shock  cannot  be  thicker  than  the  shock 
layer,  i.e,  §  A  0(A)  .  Under  these  restrictions,  the  errors  in  the  analyses 
may  be  kept  at  essentially  the  same  order  as  £  ,  or,  at  most,  LjlwL  . 

The  above  condition,  $  c.  0(A)  ■  may  be  expressed  with  the  help  of 
estimates  based  on  the  Navier-Stokes  equations  in  terms  of  the  parameter  K 
previously  introduced.  The  essential  requirements  for  the  present  formulation 
may  then  be  written  as 

6  24. f  ,  K*  >  0(0 

which  are  consistent  with  Eq.  (3.2)  of  Part  I.  Implicit  he  rein,  of  course,  is 

also  the  hypersonic  condition  t  i  .  With  some  sacrifice  in  accuracy, 

the  present  formulation  could  also  be  extended  to  cover  completely  the  fully  - 

merged-layer  regime  of  Probstein  and  Kemp.  ^  At  the  lower  Reynolds-number 

/  _ —  > 

end  of  the  latter  regime  VC5  yi  d)  ;  the  errors  in  the  present  theory  would 
become  /6~ JU>  t  instead  of  £  jlr\/  *  . 


For  the  surface  pressure  to  be  correct,  one  requires,  in  addition  to  a  unit 
surface  accommodation  coefficient  for  the  normal  component,  that  Tw /  Ta  I 


For  pointed  bodies,  the  quantity  CL  in  <  —  £.  — ° 


*  ^30  ^CO  ^ 


must  be  re- 


/“*  T0 

placed  of  course  by  a  length  representative  of  the  lateral  dimension  of  the  body. 
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Thus,  using  an  approximation  consistent  with  the  thin- shock-laye r  con¬ 
cept,  one  can  determine  the  hypersonic  viscous  flow  field  around  a  blunt  body 
independent  of  the  structure  of  the  shock-transition  zone,  even  though  the  shock 
and  shock  layer  are  comparable  in  thickness.  The  above  formulation  clearly 
indicates  that  insofar  as  the  skin  friction  and  heat-transfer  characteristics  are 
concerned,  the  shock-wave  structure  and  related  thickness-curvature  effects 
are  only  secondary  in  importance.  For  purpose  of  analysis,  the  shock-transition 
zone  may  be  treated  as  a  surface  of  discontinuity  located  at  the  outer  edge  of 
the  shock  layer. 

As  do  the  boundary-layer  equations,  the  partial  differential  equations  Eq. 
(5.22)  lead  to  an  initial-value  problem.  With  U.  ,  H  and  -f*  specified  at  any 
station  /£  ,  the  solution  can  be  continued  downstream  by  a  forward  integration 

with  respect  to  #  ,  satisfying  also  the  boundary  conditions  Eos.  (5.23)  and 

(5.24).  The  method  of  solution  to  this  system  of  equations  will  be  discussed  and 
developed  in  Section  6;  the  theory  will  be  applied  to  study  flows  around  simple 
blunt  shapes  in  Section  7. 

5.  6  Assessment  of  the  Accuracy 

Before  going  on  to  the  applications  of  Sections  6  to  8,  some  indication  of 
the  accuracy  of  the  present  approach  is  desirable  particularly  in  view  of  its  sim¬ 
plicity.  For  this  purpose,  certain  results  obtained  by  application  of  the  present 
theory  will  be  singled  out  for  discussion  here.  Two  comparisons  are  made  with 
the  results  of  previous,  more  exact  analyses.  In  one,  the  shock-layer  theory  of 
the  present  formulation  is  applied  in  the  inviscid  limit  to  describe  the  hypersonic 
flow  field  and  body  supporting  a  paraboloidal  shock;  the  results  are  then  com¬ 
pared  with  the  corresponding  "exact"  numerical  solutions  of  Van  Dyke^  and 
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53 

Swigart.  In  the  other,  the  present  theory  is  compared  with  a  numerical  solu- 

40 

tion  of  Levinsky  and  Yoshihara  for  the  stagnation  region  in  the  merged-layer 

regimes.  The  comparison  in  the  first  case  bears  on  the  general  criticism  with 

regard  to  the  poor  degree  of  approximation  of  the  Newtonian  shock-layer 
1  4  12  55 

theory,  ’  ’  ’  as  discussed  in  Section  5. 2,  The  example  chosen  is  particu¬ 

larly  meaningful  as  it  has  a  direct  implication  for  accuracy  of  the  subsequent 
examples  of  Section  7.  The  comparison  in  the  second  case  provides  a  crucial 
test  of  the  present  two-layer  model 

Paraboloidal  Shock  in  the  Inviscid  Limit  -  The  inviscid  supersonic  flow  asso¬ 
ciated  with  a  paraboloidal  shock  has  been  analyzed  numerically,  for  M  ' — *00  , 

1  53 

Y  -  1.40,  by  Van  Dyke  and  Swigart.  For  the  assessment  of  the  accuracy 

of  the  thin- shock-layer  formulation  based  on  Eqs.  (5. 22)  to  (5.  24),  the  same 

problem  has  been  analyzed  in  the  limit  of  a  vanishing  viscosity  ju. — 0  ,  using 

Eqs.  (5  .  1  3)  and  (5  .  1  5).  The  non- slip  boundary  conditions  Eq.  (5.14)  are  not 

used,  of  course.  The  results  of  the  analysis  are  presented  in  Fig.  5.  2  along 

1  .53 

with  the  corresponding  results  of  Van  Dyke  and  Swigart.  To  avoid  distraction 
from  the  basic  point  of  interest  here,  the  method  employed  to  solve  Eqs.  (5.22) 
and  (5. 24)  will  not  be  discussed  except  to  mention  that  it  was  similar  to  the 
finite  difference  scheme  described  in  Section  6. 

Figure  5.2(a)  gives  the  streamline  patterns  and  the  body  shapes  support¬ 
ing  the  paraboloidal  bow  shock.  Figure  5.2(b)  presents  the  pressure  distributions 
along  those  body  surfaces.  Unlike  the  corresponding  Newtonian  theory  (results 
not  shown),  the  present  formulation  shows  reasonably  good  agreement  with  the 
"exact"  solution.  The  body  surfaces  predicted  by  both  the  exact  method  and  the 
present  approximation  agree  very  well.  The  present  calculation  provides  a 
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somewhat  higher  pressure  at  the  surface;  the  difference  of  the  two  pressure 
distributions  varies  from  6%  at  the  stagnation  region  to  about  18%  beyond  the 
sonic  region.  Note  that  if  the  pressure  were  plotted  in  the  form  of  the  pressure 
ratio  -p>  0  .  the  agreement  would  be  even  better. 

One  notes  in  Fig.  5.2(a)  a  discrepancy  in  the  location  of  sonic  lines  in  the 
two  analyses.  This  can,  however,  be  explained  by  the  observation  that  the 
changes  in  Mach  number  along  the  streamlines  in  the  shock  layer  are  generally 
quite  slow  (the  changes  being  of  the  order  /£  according  to  the  Newtonian 
theory).  The  differences  in  Mach  number  at  the  same  point  are  actually  quite 
small,  being  typically  15%. 

The  comparison  given  above  shows  that  the  present  version  of  the  thin- 
shock-layer  analysis  is  reasonably  adequate  for  the  study  of  flow  fields  assoc¬ 
iated  with  a  paraboloidal  or  similar  shock.  The  results  have  also  confirmed 
the  importance  of  the  tangentiaL  pressure  gradient  effect  in  the  analysis  of  the 
hypersonic  shock  layer  (cf.  Section  5.2). 

Merged-Layer  Stagnation  Region  -  The  analysis  of  the  hypersonic  stagnation 

40 

region  in  the  merged-layer  regimes  by  Levinsky  and  Yoshihara  has  been 
discussed  in  Section  3  of  Part  I.  The  particular  set  of  results  of  Ref.  40 
singled  out  for  the  purpose  of  comparison  here  was  calculated  for  the  case  of 
=  10,  Y  =  5/3  and  Re  e  =  100,  corresponding  to  £  =  1/5  and  =  5. 

In  their  analysis,  the  Prandtl  number  was  assumed  to  be  3/4,  the  wall  temper¬ 
ature  taken  to  be  equal  to  the  free- stream  value,  and  the  viscosity  law  yU  oC  /y~ 
adopted.  The  results  for  the  profiles  of  temperature,  normal  velocity  and 
tangential  velocity  are  reproduced  in  Fig.  5.3  as  solid  curves.  To  distinguish 
from  their  results,  the  corresponding  analyses  based  on  the  two-layer  model 
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described  in  Section  5.  4  and  5.  5  are  represented  in  the  same  figure  as 
circles  and  triangles.  Considering  the  value  T"  =  1.667,  the  applicability 
of  the  thin- shock-layer  theory  to  the  present  problem  should  be  regarded  as 
being  critical. 

The  present  analysis  gives  a  shock  interface  at  'ty/0'  &  0.10.  The 

solution  for  /(L  0.  10  is  obtained  by  numerical  integration  of  the  shock- 

layer  equations  of  Eq.  (5.22)  to  (5.24)  in  the  stagnation  region,  and  the  solution 
for  0.  10  is  obtained  from  the  transition- zone  solution  of  Eqs.  (5.  15)  to 

are  provided  by  the  shock-layer 


(5.  18).  The  value  of  Uz  ,  Hz  and 


^'‘'Ab¬ 
solution.  In  addition  to  the  profiles  of  the  normal  velocity  /\r  ,  the  temperature 
T  and  the  tangential-velocity  gradient  U  which  are  reproduced  from  Ref.  40, 
the  corresponding  density  distribution  p/ f3^  in  the  shock-transition  zone  is  also 
calculated  and  presented  in  Fig.  5.  3  for  comparison. 

As  is  evident  from  the  figure,  the  analysis  of  Levinsky  and  Yoshihara 
and  the  considerably  simplified  treatment  based  on  the  present  formulation 
agree  in  all  essential  details  both  in  the  shock  layer  and  in  the  shock-transition 
zone.  The  small  difference  between  the  two  analyses  is  somewhat  surprising, 
when  one  notes  that  the  specific-heat  ratio  1.667  is  rather  far  from  unity.  One 
notes,  however,  that  the  profiles  of  normal  velocity  and  temperature  given  by 
the  present  analysis  are  discontinuous  at  the  shock  interface  0.10. 

These  discontinuities  are  unavoidable  in  the  present  formulation  which  matches 
the  derivatives  of  ^  at  the  interface.  However,  they  are  rather  small,  being, 
of  order  £  as  compared  to  unity,  and  are  indeed  consistent  with  the  assumption 
of  the  theory  (see  discussion  of  the  preceding  section). 

It  may  be  noted  that  if  one  were  to  interpret  the  length  "a,"  in  the  shock- 
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layer  solution  as  the  shock  radius  instead  of  the  body  nose  radius,  as  was  done 
in  the  preceding  paraboloidal- shock  problem,  the  abscisa  of  the  circles  and 
triangles  in  Fig.  5,  3  would  increase  by  10%  and  result  in  a  better  agreement 
with  the  solid  curves.  This  interpretation  is  permissable  within  the  present 
thin-layer  approximation  and  also  implies  that  the  shock  and  the  body  are  con¬ 
centric.  The  latter  assumption  is  in  fact  implicit  in  the  analysis  of  Ref.  40. 

It  must  be  pointed  out  that  the  reasonably  good  agreement  revealed  in  Fig. 

5.  3  is  by  no  means  an  indication  of  the  accuracy  of  present  formulation  in  the 
absolute  sense,  since  the  analysis  of  Levinsky  and  Yoshihara  also  involves  cer¬ 
tain  approximations.  The  small  differences  do  demonstrate  convincingly,  how¬ 
ever,  the  validity  and  adequacy  of  the  two-layer  concept  underlying  the  present 
formulation.  The  result  of  the  comparison  implies,  of  course,  that  the  most 
important  aspect  of  the  flow  phenomena  affecting  the  aerodynamics  and  heat 
transfer  in  the  merged-layer  regimes  is  not  the  direct  effects  related  to  the 
shock-thickness  and  the  shock-wave  structure  but  rather  the  slip-like  transport 
effects  immediately  behind  the  shock. 

6.  METHOD  FOR  SOLVING  THE  VISCOUS  SHOCK  LAYER  EQUATIONS 

With  the  basic  problem  formulated  and  the  adequacy  of  the  flow  model 
assessed,  the  thin- shock-layer  theory  will  now  be  applied  to  study  hyper  sonic 
viscous  flows  over  wedge,  paraboloid  and  other  simple  blunt  shapes.  To  do  so, 
we  must  first  develop  a  method  which  is  capable  of  obtaining  nonsimilar  sol¬ 
utions  to  the  partial  differential  equations  Eqs.  (5.22)  to  (5.  24)  governing  the 
shock  layer.  The  present  solution  to  the  dissimilar  problem  requires  the 
combined  use  of  two  methods;  a  forward -integration  method  based  on  the  finite- 
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difference  approximation,  and  the  solution  by  scries  in  the  vicinity  of  the 
leading  edge  or  the  stagnation  point.  The  iar.ier  solution  is  required  tc  pro¬ 
duce  valid  initial  data  for  calculation  by  the  finite -difference  method. 

6.1  Discussion  of  General  Treatment 

For  the  analysis  of  the  shock  layer,  it  is  convenient  to  work  with  the 
dimensionless  quantities  ,  Q  and  -p_. 


U  = 


CO 


©  = 


h  -  H 


+/ 


Hco~Hw 


P 


fvo  'JSO 


(6. 


Introduce  the  variables 

4  ^  a. 

Note  that  is  the  distance  measured  along  the  axis  of  symmetry  (cf  Fig. 

5.1)  related  to  of  through  (cL  oc)  +  (d  Z  )  >  and  that  u  and  © 

behave  like  ~/y  near  the  body  surface  and  therefore  will  be  regular  in  £ 

In  terms  of  the  new  variables  of  Eq.  (6.1)  and  (6.2),  the  system  of  differential 
equations  Eqs.  (5.  22}  is  transformed  to 
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The  subscripts’ A  and  4  denote  partial  derivative  with  respect  to  A  and 
£  .  As  before,  the  parameter  l(  is  defined  as 

W  1  —  r  '°<x>  ’  <#  ^  c  O  i  y^o  Ttf  ^ 

K  =  £  ~~7*7~  XI '  -  ^ '(XT' ' T7 

where  T  and  are  the  reference  temperature  and  the  corresponding 

viscosity  respectively.  If  the  linear  representation  =  is  adopted, 

the  factor  l\j  in  the  above  equations  becomes  unity.  The  boundary  conditions 
at  the  surface  2^=0,  Eq.  (5.23),  becomes 

£  =  &  -Q  (6.4) 

The  shock  relations  at  the  outer  edge  of  the  shock  layer,  Eq.  (5.24),  become 
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5  H) 

These  relations  are  now  to  be  satisfied  at  £  =  1  .  The  relation  between  the 
physical  distance  from  the  body  surface  and  £  can  be  determined  from 


JL  -(  2  \ 


(Y< 


£  2  "  y  <r*t,/T 


j  Tv 

IT 


T» 


a* 

— =—  —  CO<*  p 
U  ! 
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& 


(6. 6) 


The  local  skin  friction  and  the  local  heat-transfer  rate  at  the  surface  can  be 

calculated  from  p  i  (A  -  and  0  _  at  £  =0: 

£  ^ 


'  J  <  v,r  | 

(/*  ^rL 
e r  s  -V 

2  ,Pco  '  4o 


O'  £*y~  |  /  . 

V*/  ' 


c«  5 


§T-\ 

(j  ■ —  i 

A _ 


(6.7) 


_  I  f-  J 
2  <■“ 


c ■p'  >  & 


p  U  'H  ~  Hj  "  2R-<*  "v_  r/  r  ? 

Aside  from  the  dependence  of  the  solution  on  the  shock  or  body  shape 
through  ft  and  CT  ,  the  system  Eqs.  (6.  3)  to  (6.  5)  is  governed  by  the  par¬ 
ameters:  «2  ,  £  ,7 v /T0  and  pr-  .  Clearly,  <"  *  is  the  most  important 

parameter  and  controls  the  transport  effects  in  the  shock  layer.  Following 
the  classification  for  the  stagnation  region,  Eq,  (3.2)  of  Part  I,  one  may  also 
divide  the  general  shock-layer  problem  into  two  regimes,  in  one,  defined  by 
0(D  £  £K2  4  oo  ,  the  inviscia  sh  ock  r  elations  hold,  i.e,  U — >© — >-p> — >j 

at  £,  =  1  .  In  the  other,  where  0\6)  —  £■  K  —0(i)  .  the  speed  and  enthalpy 
changes  across  the  shock-transition  zone  are  important. 

As  noted  before,  the  thin- shock-layer  approximation  has  reduced  the 
governing  differential  equations  to  the  parabolic  type.  This  makes  forward 
integration  in  the  downstream  direction  possible.  One  may  note,  however, 
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that  the  great  simplification  found  in  the  analysis  of  the  stagnation  region 
cannot  be  found  in  the  general  case.  Specifically,  an  analytical  treatment 
in  the  merged-layer  regimes  (Regime  II)  is  not  generally  possible.  Although 
the  differential  equations  themselves  admit  a  separation  of  variables  (after 
neglecting  the  tangential-pressure-gradient  term),  a  self-similar  solution 
does  not  exist  owing  to  the  outer  boundary  condition  Eq.  (6.  5). 

For  analysis  in  the  higher  Reynolds  number  regime  (Regime  I),  Eqs. 

(6.  3)  to  (6. 5)  may  best  be  handled  by  transformation  to  the  coordinate  system 
(  familiar  inboundary-layer  theory.  However,  vorticity-interaction 

theory  based  on  the  boundary-layer  approximation  can  no  longer  be  extended, 
as  is  permissible  for  the  stagnation  region,  to  treat  the  non-linear  vorticity 
effect  in  the  viscous-layer  regime.  This  is  because  the  solution  in  the  more 
general  case  cannot  be  simply  described  in  terms  of  the  asymptotic  property 
of  the  boundary-layer  solution,  as  in  the  case  of  the  stagnation  region  (or,  at 
least,  the  possibility  of  such  a  simple  description  is  not  obvious).  Of  course, 
for  small  departures  from  the  boundary-layer  limit,  the  second-order 
boundary-layer  theory  is  valid.  In  fact,  Hayes'  original  matching  condition, 
Eq.  (1.1) 

f*  =  j  +  zn  e  H — ~ 

is  consistent  with  the  thin- shock-laye  r  approach  and  should  be  generally 
applicable.  A  treatment  of  the  vorticity-interaction  theory  based  on  Eq. 


Presumably,  the  alternative  condition  Eq.  (1  .  3)  ■P„„(oo)  =  Ci  may  also 

apply.  To  the  writer 1  s  knowledge  ,  the  equivalence  oi  Eqs.  ( 1  .  1 )  and  ( 1  .  3) 
has  not  been  proved  in  the  general  case,  however. 
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(1.1)  (which  also  requires  the  use  of  a  numerical  method  to  handle  the  non¬ 
similar  solution),  though  not  carried  out  here,  may  be  considered  as  com¬ 
plementary  to  the  method  presently  to  be  described,  whose  range  of  validity 
covers  the  merged-layer  as  well  as  the  viscous-layer  regimes 

Before  passing  to  the  discussion  of  the  two  methods  employed  in  the 

4  29  61 

present  procedure,  the  possibility  of  using  momentum -integral  methods  ’  ’ 

6>  Z  63 

as  well  as  other  methods  familiar  from  the  boundary -layer  theories,  ’ 
should  not  be  completely  overlooked.  The  main  objection  to  the  use  of 
momentum-integral  methods  which  also  involve  stepwise  numerical  integration 
for  the  present  study  is  its  limited  range  of  applicability.  This  shortcoming  is 
related  to  the  fact  that  the  velocity  and  temperature  profiles  of  the  shock  layer 
in  the  low  and  high  Reynolds  number  ranges  considered  are  radically  different, 
and  do  not  permit  the  simple  description  characteristic  of  most  integral  methods. 
6.  2  Series  Expansion  Near  Axis  of  Symmetry 

For  flow  regions  not  far  removed  from  the  axis  of  symmetry,  solution 
by  a  series  in  ascending  power  of  the  distance  from  the  axis  is  possible.  Thus 
for  the  stagnation  region,  one  may  assume  a  development  in  ( ' %/ 1)  for  u  , 

@  and  -fv  .  In  terms  of  the  variables  A  and  ,  one  has,  for  the  stagna¬ 
tion  region  of  a  smooth  blunt  body  (considering  only  positive  X  ). 


U  -  tL  +  A  M.  +  \2  a  .  + .  (6-8) 

with  similar  developments  for  0  and  -p;  ,  The  leading  terms  Ua  ,  Q0 

and  ,  together  with  -p>  ,  determine  the  solution  in  the  stagnation  region 
for  which  Eqs.  (6.  3)  to  (6.5)  reduce  to  the  following  system  of  ordinary 
differential  equations 
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'P'l  =0 

f'',  =(2-^)aeC  y 


\  (0)  =  0/0)  =  0 
fi'o(l)  -  I  >  -p/,(l)~ 0 

ujl)  =  1  --J^r  i-ItJOSjO) 

8,M  =  i-5^Y  N  Mue.'fO 

where  the  prime  (  )'  denotes  differentiation  with  respect  to  £  .  The 

appearance  of  in  the  first  equations  signifies  the  tangential-pressure- 
gradient  effect.  Eqs.  (6.9)  to  (6,  11)  constitute  a  two-point  boundary-value 
problem.  This  analysis  of  the  stagnation  region  has  been  discussed  in  Section 
3  of  Part  I.  More  detailed  results  are  given  in  Refs.  42-44. 

By  collecting  terms  of  the  next  power  of  \  in  the  Eqs.  (6.  3)  to  (6.  5), 
one  obtains  a  system  of  ordinary  differential  equations  governing  W,  ,  ®  , 

and  A'  These  equations,  given  in  the  Appendix,  lead  to  a  two-point 
boundary-value  problem  involving  a  system  of  linear  ordinary  differential 
equations.  By  virtue  of  their  linearity,  however,  solution  to  this  two-point 
boundary- value  problem  can  be  obtained  simply  by  superposition  of  the  par¬ 
ticular  solutions.  This  considerably  simplifies  the  task  of  numerical  analyses. 

The  systems  of  equations  governing  the  coefficients  for  the  higher  powers 


\ 


(6.10) 


(6.11) 
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°f  A  >  U  y  ’  an<^  “fX/3  '  etc'  are  Slmlfar>  but  more  complex.  In  the 

following  study  of  the  problems  of  paraboloidal  and  hyperboloidal  shapes, 
only  the  coefficients  U.  ,  ©/  and  -p/^  have  been  determined.  The  range  of 
validity  of  the  solution  based  on  the  two-term  representation,  Un+\Hi  ,  etc., 
is  of  course  rather  limited,  but  the  roles  of  the  coefficients  >  ©(  and  A 
are  nevertheless  crucial  for  the  numerical  solution  of  the  downstream  flow 
field,  as  will  be  discussed  subsequently. 

The  series  solution  of  the  type  of  Eq.  (6.  8)  is  applicable  only  for  a  body 
with  finite  nose  radius.  F or  a  pointed  body,  the  appropriate  series  solutions 
take  the  form 

u  —  u.o  A  +  u  t  A  +  uz  A  +  •  •  •  • 


0  =  &  X'/S  +  0,4  -f  ©s  X'*  + 
_ 

•f'  -  /  -h  A  }  +  A  f 


(6.12) 


The  leading  terms  for  U  and  0  are  simply  proportional  to  £  .  In  this 

case,  the  successive  coefficients  U ,  ,  U  2  ,  ,  etc.  can  also  be  determined 

in  explicit  analytical  form.  These  developments  will  be  used  in  the  study  of 
wedge  and  cone  flows  in  Section  7.1. 

6.  3  Finite -Difference  Method 

Numerical  methods  for  solving  the  parabolic  and  hyperbolic  partial 
differential  equations  by  finite-difference  approximations  have  been  quite 
extensively  studied  in  recent  years  (for  example,  see  Refs  64  and  65).  The 
methods  have  been  applied  to  viscous -flow  problems  within  the  context  of  the 
boundary-layer  theory  by  a  number  of  investigators.  ^  ^ 

Singularities  in  the  Difference  Equations  -  In  the  finite-difference 
approximation,  the  derivatives  are  replaced  by  the  quotients  of  the  differences, 
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and  the  partial  differential  equations  are  satisfied  approximately  at  a  finite 
number  of  grid  or  lattice  points.  For  an  initial- value  problem  such  as  the 
present  one,  a  forward  integration  is  possible;  the  solution  to  the  difference 
equations  may  then  be  determined  in  terms  of  values  obtained  at  preceding 
stations  or  grid  lines.  That  is,  the  difference  approximation  is  used  to 
evaluate  dQjdX  in  terms  of  quantities  given  on  the  right-hand  side 

of  Eq.  (b.  3).  The  quantities  U  and  0  are  then  obtained  by  integrating  with 
respect  to  \  .  Strictly  speaking,  this  procedure  is  not  applicable  to  Eq. 

(6.  3)  at  X  -  0  and  ^  -  0.  because  some  of  the  terms  on  the  right-hand  side 
of  the  first  and  second  of  Eqs.  (6.  3)  are  proportional  to 


\*04,j6)  ~  \C 

Hence  in  the  limit  of  either  X — 0  or  £,—> 0  ,  the  values  of  Su/$\  ,  etc. 

become  indeterminate,  or  even  infinite  (in  the  case  of  pointed  bodies) .  This 

leads  to  two  important  observations.  First,  in  order  to  achieve  forward 

integration  correctly  in  the  vicinity  of  the  stagnation  point  X  =  0,  it  is  most 

essential  to  provide  the  coefficients  ,  etc.  in  the  series  solution  of  Eq. 

(6.  8),  which  are  in  fact  the  correct  expressions  for  ,  etc.  at  X  =  0. 

Starting  the  integration  procedure  merely  with  the  knowledge  of  U0  ,  etc.  , 

66  67 

as  has  been  proposed  in  some  previous  work,  ’  cannot  yield  the  correct 
solution.  Second,  in  view  of  the  factor  f/AC  ,  errors  resulting  from  the 
difference  approximations  will  be  greatly  amplified  in  the  vicinity  of  the 
"leading  edge"  X  =  0  and  at  the  body  surface  C  -  0.  As  far  as  the  method 
of  difference  is  concerned,  the  lines  \  =  0  and  =0  may  therefore  be 

regarded  as  singularities  of  the  equations.  It  is  clear  that  application  of  the 
difference  equations  must  exclude  these  lines. 


Ll 


I 
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The  manner  in  which  the  difference  quotients  are  formed  and  in  which 
the  equations  are  approximated  gives  rise  to  a  variety  of  difference  methods. 

The  best  known  are  the  standard  explicit  and  implicit  scheme  and,  in  addition, 
that  introduced  by  Dufort  and  Frankel.  Near  the  singularities,  the 

conditions  for  stability  and  convergence,  together  with  the  small  truncation- 
error  requirement,  are  so  stringent  that  the  explicit  schemes  become  hardly 
practical  for  the  present  application.  The  implicit  scheme,  which  is  sub¬ 
ject  to  none  of  these  restrictions,  not  only  is  superior  as  has  been  demonstrated 

6  7 

by  Flugge-Lotz  and  Blottner  for  the  boundary-layer  analyses,  but  is 

necessary  in  the  present  analysis.  This  is  quite  clearly  brought  out  by  the 

71 

conclusion  of  a  recent  study  of  A.  L.  Chang  and  the  writer,  in  which  the 
relative  merit  and  efficiency  of  all  three  schemes  are  examined  for  the  case 
of  a  cone  (the  result  of  which  is  to  be  presented  in  Section  7  ). 

In  view  of  the  critical  nature  of  the  singularities  in  amplifying  the  trun¬ 
cation  errors,  the  regions  near  =0  and  £  =0  will  both  be  excluded.  The 

system  of  difference  equations  considered  applies  therefore  only  in  the  domain 


*  1 

;.4c  i 


The  boundary  condition  originally  at  U  (0)  -  0  (0)  =  0  can  be  simply  trans¬ 
formed  to  the  new  boundary  =  ~C,  by  the  Taylor  theorem 

-  (6.14) 


The  stability  for  the  standard  explicit  scheme  would  require  dA  -  0  (KfoO’JL 
the  convergence  of  the  Dufort-Frankel  scheme  would  r  equir  e  (AA.)Z  • 

Whereas  small  trunction  errors  would  require  $  «i 
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where  U  ,  U  ^  and  are  evaluated  at  £  =  •  If  one  applies  the 

first-order  correction  in  the  Taylor  series,  as  is  done  below,  the  error  in¬ 
curred  by  shifting  the  boundary  from  £  =  0  to  £  ■=  will  be  of  the  order 

2,  — 

£  .  In  the  following,  the  initial  data  for  LL  ,  etc.  will  be  prescribed  at 

A0  by  the  two-term  development  >  etc.  The  exclusion  of  A  -  ^  —  0 

in  this  manner  will  therefore  introduce  an  error  of  the  order  (In 

most  calculations  performed  Xo  and  are  taken  as  0.  10.  )  With  the 
exclusion  of  the  region  close  to  A  =  £  -  0  ,  the  truncation  errors  in  the  dif¬ 


ference  equations  are  still  very  large;  solution  with  acceptable  accuracy  will 
require  increments  A  X  and  A  £>  considerably  smaller  than  and  £  . 

(In  subsequent  applications,  A  A  and  A  ^  are  taken  in  most  cases  as  1/100.) 

Difference  Quotients  and  Difference  Equations  -  In  the  implicit  scheme 
(see  the  sketch  on  Fig.  6.  1)  the  unknown  values  of  U  and  ©  at  point  e 
along  a  new  grid  line  (column  m)  are  to  be  solved  simultaneously  with  values 
at  points  d,  f,  etc.  along  the  same  grid  line  using  known  data  from  the  pre¬ 
ceding  stations.  The  difference  quotients  for  U  derivatives  which  define 
the  present  scheme,  are  the  first  terms  in  the  following  development 

taTT.  - -  f  -T&A)  . 


2(AA) 


Uf~Ud.  __L/a  £  +  , 

2(A£) 


_ 


(6  14) 


To  be  more  precise,  the  following  formulation  requires  initial  data  along 

A=A„  as  well  as  along  A  =  A"AA  • 
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The  difference  quotients  for  Q  derivatives  arc  similar.  One  must  note 
that  the  difference  approximations  for  3u/d\  have  been  carried  out  to  within 
an  error  of  order  (A  A)2  instead  of  ^A  X)  ■  With  these  quotients,  the  partial 
differential  equations  Eqs.  (6.  3)  are  approximately  satisfied  at  the  grid  point 
e  .  The  difference  equations  can  be  linearized  without  sacrificing  accuracy 


by  expressing  the  values  of  the  coefficients  of 


3  u.  3*  u 


j-  ,  etc .  in  the  dif- 


3C  ’ 

ferential  equations  at  point  ^  in  terms  of  their  corresponding  values  at  b 
and  of  the  preceding  station  (refer  to  Fig.  6.1).  This  is  accomplished 
through  the  relation 


Xe  =  (2Xb-XJ- 


(2U) 


<X) 


e 


(6.  15) 


where  X  may  be  any  of  the  dependent  variables.  To  eliminate  the  pressure 
as  an  unknown  in  the  difference  equation,  one  can  replace  in  terms  of 

values  of  at  the  preceding  three  stations  b  ,  and  Jv  (see  sketch) 


through 


££)  =  2 
3*  4 


3A 


_  ■*  3  *<.  +  o  feu) 

2.  (AX) 


e/AAA 
2. 


+  -  . 


(6. 16) 


The  pressure  -ps  can  always  be  obtained  after  u  is  determined  by  integrating 
the  third  of  Eq.  (6.  3)  with  respect  to  ^  .  The  error  resulting  from  the 

integration  by  Simpson's  rule  is  of  the  order  Application  of  Eqs.  (6.  14) 

to  (6.  16)  to  the  system  of  differential  equations  Eqs.  (6.  3)  gives  rise  to  twice 
as  many  algebraic  equations  as  the  number  of  points  between  -  C  and 
along  a  grid  line. 

To  complete  the  formulation,  one  must  convert  the  boundary  condition 
at  £  =1,  Eq.  (6.5),  and  the  boundary  condition  at  £  =  ,  Eq.  (6.13),  to 


AF  - 1  285-A- 1 0 


56 


two  linear  algebraic  relations  among  the  unknown  values  of  U  and  0  without 
introducing  new  unknowns  corresponding  to  points  beyond  £^4.  £  <  /  •  This 

can  be  accomplished  by  the  typical  relation 

(I t)f=  ZST7[3  V4a*+a»]  +0(A?f  <£"17) 

and  relations  similar  to  Eq.  (6.  15).  The  difference  equations  at  all  grid 
points  along  the  same  grid,  together  with  the  boundary  conditions,  form  a 
linear  algebraic  system  sufficient  for  determining  W  and  @  at  all  grid 
points  between  and  ^  =  1  .  The  complete  system  of  equations  is  given 

in  the  Appendix.  It  sufficies  to  say  that  the  matrix  of  this  system  is  the  tri¬ 
diagonal  type,  and  the  solution  can  be  obtained  by  following  the  standard 

6  5 

procedure,  e.g.  Gauss'  elimination  method. 

Truncation  Error  -  One  may  now  examine  more  closely  the  truncation 
errors  in  the  difference  equations.  Only  the  truncation  errors  associated 
with  the  first  of  Eq.  (6.  3)  need  be  considered.  A  similar  conclusion  applies 
to  the  other  equations.  There  are  two  terms  in  the  remainders  of  the  dif¬ 
ference  equations  which  are  most  critically  affected  by  the  singularity 
A  =0  and  C,  =  0:  one  associated  with  the  tangential  pressure  gradient  and 
the  other  with  the  tangential  stress  term.  With  the  knowledge  of  the  remainders 
given  in  Eqs.  (6.  14)  to  (6.16),  and  of  the  fact  that  near  5  -  Q>  A 

and  U  £  ,  the  fractional  error  in  u.  caused  by  these  two  terms  may  be 


estimated  as 

.  Tw  (AA)Z 


A?2 


u  + 

55^ 


e3L‘s 

K z  « 


The  factors 


are  not  large  numbers,  especially 
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at  K  2  -  0(l)  >  because  U.  is  essentially  proportional  to  C,  for  small 
X  (refer  to  the  discussion  in  Section  3,  Part.  1)  In  the  merged-laver 
regimes,  at  least,  the  fractional  error  in  U  cannot  exceed 

6U.  .  Tw  (A  A)*  . 

~  — 7-7—  +  - =; — F -  > 

since 

A-  ^  A  3  .»  &  > 

Taking  XQ  ^  an<a  A  X  A  t,  .«o  ■y——  .  and  £  ^  -  ,  the  fractional 

error  S  U  /  U.  is  seen  to  remain  at  the  level  of  one  percent, 

7.  STUDY  OF  FLOY/S  OVER  SIMPLE  NONSLENDER  SHAPES 

The  method  of  solution  formulated  in  the  preceding  section  will  now 

be  applied  to  study  hypersonic  flow  around  wedge,  cone,  paraboloid  and 

hyperboloid  shapes  at  low  Reynolds  number  Only  solutions  to  the  shock- 

layer  region  will  be  discussed  below.  The  flow  structure  in  the  shock- 

transition  zone,  in  each  case,  can  be  determined  by  simply  integrating 

the  first-order  ordinary  differential  equations,  Eq,  (5.13),  after  the  shock- 

layer  solution  is  obtained. 

7.1  Nonslender  Cone  and  Wedge 

As  examples  of  flows  around  pointed  bodies,  the  shock  layers  on 

cones  and  wedges  will  first  be  studied.  The  analysis  may  provide  a  basis 

for  assessing  the  low-density  aerodynamic  and  heat-transfer  characteristics 

of  flat  surface  at  high  incidence,  which  could  be  meaningful  to  the  current 

72-74 

lifting  re-entry  studies.  From  the  viewpoint  of  present  approach 

based  on  the  two-layer  concept,  the  cases  considered  are  interesting  in 
that  the  curvature  of  shock  above  the  wedge,  and  cone  surfaces  is  generally 
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small,  and,  therefore,  the  formulation  developed  in  Section  5.  5,  which  has 
neglected  the  shock  curvature  effect,  could  possess  a  regime  of  validity 
even  beyond  the  limit 


^  o(s) 


On  the  other  hand,  the  solution  will  break  down  in  the  vicinity  of  the  leading 
edge  because  of  its  singularity  at  the  apex.  This  will  be  discussed  later. 

One  may  also  note  that  the  problem  of  a  flat-plate  at  zero  incidence 
has  been  the  subject  of  both  theoretical  and  experimental  studies  of  recent 
years.^^’  ^  ^  Generally  speaking,  the  problem  is  complicated  by  the 
presence  of  both  the  boundary-layer  displacement  and  surface- slip  effects. 
However,  in  the  present  problem  which  considers  thin  shock  layers  over 
highly  cooled,  nonslender  cones  and  wedges,  these  two  effects  can  be 
regarded  as  unimportant  compared  to  the  slip-like  velocity  and  enthalpy 
changes  across  the  shock. 

The  problem  under  the  present  framework  is  also  attractive  in  that 

for  the  case  of  a  linear  viscosity  law  the  analysis  is  reduced  to  the  solution 

of  a  single  parabolic  differential  equation  governing  the  velocity.  It  may 

therefore  serve  as  a  testing  ground  for  the  various  methods  and  schemes 

71 

of  solving  the  shock-layer  equations.  In  this  case,  the  first  of  Eqs. 

-1/ 


(6.  3)  can  be  alternatively  written  for  YJ=U/Z  as, 

S-SsT 


31V  +(,+  jf  yzw 


5T 


H  3T' 


(7.1) 


The  boundary  conditions,  Eqs.  (6.4)  and  (6.5),  are 
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W  =  0 


yzw  =  i-1?- 

'fj  <3  if 


yr  =  o) 


(r  =  /) 


Note  Eq.  (7.1)  contains  none  of  parameters  5  ,  £  ,  pr  and  Tw /T0  • 

The  effects  of  gas  rarefaction  and  surface  incidence,  as  well  as  the  variation 
of  the  reference  temperature  with  respect  to  /£  ,  are  all  contained  in,  and 


therefore  can  all  be  correlated  by,  the  variable 

The  appearance  of  the  factor  in  v  signifies  the  influence  of  the 

oblique  -  shock  angle  on  the  local  Reynolds  number.  This  is  particularly 
evident,  if  one  considers  a  linear  viscosity-temperature  relation  and  high 
Reynolds  number,  in  that  case  the  local  velocity  varies  like  /Cc >■)'/$,  and  the 

3>  — ■ 

viscosity  like  /Uhl)  .  To  relate  the  variable  ^  to  the  Reynolds  number 
based  on  the  free -stream  condition,  say  ,£<e>.  ,  one  may  use  the  reference 


(7.3) 


temperature 

Jk  =  _l  f 

T0  <7(0  j 


7"  /  T •  \  u(!)  * 

■  O  1  *  \  I  3  /  ~  3 


After  the  velocity  field  is  determined,  the  other  flow  quantities  can  in 
turn  be  found.  For  a  unit  Prandtl  number,  the  total  enthalpy  may  be  obtained 
simply  from  the  Crocco  relation 

e  -  U  =  /2W  (7 

As  observed  in  Section  6.  2,  the  solution  of  U.  may  be  developed  as 
a  series  in  the  form 


_  _  _  '/a  _  _  —  _  JA 

tA  =  u0  <£  +  U,  -fc  +  Uz  i.  + 


(7.6) 
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of  which  the  coefficients  can  be  analytically  determined,  for  -  0  and  1  , 

respectively,  as 

S0=/S"f,  C.  (V=  o) 

«.*£.  «,  =  "2  =  74?o( W7  5-2flcn^£T}tc.6M) 

The  appearance  of  the  argument  ^  ^clearly  indicates  that  the  solution  to 
the  shock-layer  equations  in  this  case  is  singular  at  the  apex  ^  =  0.  As 

in  the  classical  boundary-layer  theory,  the  immediate  vicinity  of  the  leading 
edge  of  a  cone  or  flat  plate  has  to  be  excluded  from  the  present  analysis. 

With  the  initial  data  provided  by  the  series  Eq.  (7.6),  the  numerical 
solutions  to  Eq  (7.1)  are  obtained  for  the  wedge  and  cone.  The  procedure 
used  follows  essentially  that  discussed  in  the  preceding  section,  which  is  based 
on  an  implicit  difference  scheme.  The  advantage  of  using  the  variables  W  ,  •$£■ 
and  T  rather  than  u.  ,  &  and  £,  is  that  in  these  new  variables  the 
truncation  errors  of  the  difference  scheme  for  this  particular  problem  become 
less  critical  with  respect  to  the  singularity  at  the  surface  mentioned  previously. 
Consequently,  the  increment  AV  used  is  not  required  to  be  as  small  as  A.Z,  . 
Strictly  speaking,  the  use  of  VV  and  l/T  requires  certain  corrections  to  the 
difference  quotients  shown  in  Eq.  (6.  14),  because  W  is  not  strictly  regular, 
but  behaves  like 

W  -  a  f  +  bfS/%  o(f  4) 

near  the  surface  where  -  0.  Examination  shows,  however,  that  the 

corrections  are  numerically  small  and  not  necessary,  at  least  for  the  present 


The  singularity  is  more  critical  when  -0-  this  more  general  case, 

W  =a.Y  +  b(?fZ+0‘ 
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problem . 


43 

The  results  of  skin  friction  and  heat-transfer  rate  on  cones  and  on 
wedges  for  Pf  -  1  are  presented  in  Fig.  7.  1  as 
As  shown,  the  values  of  Cw  and  approach,  respectively,  and 

^iyru2jS  as  /£, - ► 0  ,  as  required  by  the  free -molecular -flow  theory  for  unit 

accommodation  coefficients.  Fig.  7.  1  indicates  that,  as  far  as  the  heat- 
transfer  rate  and  skin  friction  on  cones  and  wedges  are  concerned,  appreciable 
departure  from  the  boundary-layer  theory  does  not  occur  until  JC  c «  /  ,  For 
4  |  ,  the  analytical  development  (6-term  expansion  for  =  0,  4-term  for 

y  =  1 )  appears  to  b  e  adequate.  The  calculation  performed  begins  at  y  =0.01 
and  uses  A.'W  =1/10  and  A  =  1/1000  for  .0/  —  •p  —‘2.  and  Ay  =1/100 
for  .2  A  —  1 0  The  total  computation  time  (including  print-out)  on  an 
IBM  704  digital  computer  is  approximately  2  minutes. 

The  corresponding  velocity  profile  tends  to  a  linear  distribution  in  5 
in  the  limit  of  — >0  ,  as  is  also  revealed  by  the  coefficients  Ua  in  the 

development  of  Eq.  (7.  5).  As  /#  increases,  the  velocity  gradient  increases 
near  the  surface,  tending  to  form  a  boundary  layer,  and  decreases  at  the  outer 
edge,  thus  reducing  the  velocity  change  across  the  shock.  The  value  of  u. 
immediately  behind  the  "shock"  is  also  given  in  Fig.  7.1.  It  is  rather  inter¬ 
esting  to  observe  in  this  respect  that,  even  at  as  low  as  3,  the  flow  field 
can  already  be  represented  very  closely  by  the  classical  boundary-layer 
limit,  in  spite  of  the  fact  that  the  velocity  at  the  outer  edge  is  still  quite  far 
from  the  inviscid  value  and  that  the  viscous  layer  is  still  a  major  fraction 
of  the  shock  layer.  To  bring  out  clearly  this  observation,  the  velocity  pro¬ 
files  at  various  stations  are  correlated  in  terms  of  the  similarity  variable 


vs.  y. 
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(7.7) 


The  good  correlation  and  the  excellent  agreement  with  the  Blasius  profile,  as 
revealed  in  Fig.  7.  2,  also  provide  a  check  on  the  accuracy  of  the  method. 

Finally,  from  the  local  skin-friction  of  the  wedge  surface,  one  can 
determine  the  lift/drag  ratio  for  a  flat  plate  at  incidence  in  the  low  Reynolds 
number  regime.  This  ratio  corresponding  to  a  46°  incidence  angle,  which 
is  calculated  on  the  basis  of  one  surface  on  the  compression  side,  is  in¬ 
cluded  at  the  lower  part  of  Fig.  7.1  as  a  function  of  .  Though  not  pre¬ 
sented  in  the  figure,  the  L/D  ratio  at  /S  =  30°  and  60°  does  not  differ  a 
great  deal  from  that  for  -  45°  .  The  L/D  ratio  at  =  1 . 0  is  about  0.35, 

becomes  lower  for  /£  <.  /  ,  and  tends  to  zero  as  <£ - > 0  ■  Once  again,  the 

present  analysis  yields  the  correct  free-molecule  limit  for  unit  accommoda¬ 
tion  coefficients.  It  may  be  noted  that  using  boundary-layer  theory  would 
lead  to  a  negative  L/D  at  low  Reynolds  number. 

As  a  check  on  the  internal  consistency  of  the  analysis,  one  can  examine 
the  solutions  to  determine  whether  or  not  the  shock  layer  remains  thin  in 
comparison  with  the  lateral  dimension  of  the  body.  Examination  of  the  solu¬ 
tion  obtained  shows  that  the  particular  assumption  of  equating  the  shock 
slope  with  the  body  slope  may  not  be  strictly  satisfactory  for  a  flat  plate 
at  less  than  25°  incidence  in  the  range  of  I  .  (The  case  of  a  cone  is 

less  critical.)  An  interesting  question  arises,  in  this  respect,  that  is, 
whether  and  when  the  viscous  shock  layer  will  also  detach  from  the  surface, 
as  it  does  in  the  inviscid  limit,  when  the  wedge  or  cone  angle  becomes  too 
high. 
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7. 2  Paraboloid  and  Hyperboloid 

As  examples  of  flows  around  bodies  of  revolution  having  finite  nose 
radii,  the  shock  layers  around  a  paraboloid  and  a  hyperboloid  will  now  be 
studied,  applying  the  full  method  developed  in  Section  6. 

The  principal  reason  for  selecting  this  class  of  bodies  is  that  the 
4  19 

zero-pressure  point  ’  does  not  occur  on  these  bodies,  and  the  shock- 
layer  theory  should  be  valid.  The  comparison  provided  at  the  end  of 
Section  5  for  the  paraboloidal  shock  in  the  inviscid  limit  has  further  indi¬ 
cated  that  the  present  version  of  the  thin- shock-layer  approximation  may 
be  used  to  describe  the  flow  field  with  fair  accuracy.  It  is  reasonable  to 
expect  that  the  approximation  will  remain  adequate  in  the  corresponding 
viscous  problem.  The  degree  of  approximation  in  the  case  of  the  hyper- 
boloidal  shape  may  presumably  be  better,  because  of  its  smaller  surface 
curvature.  Here,  the  distinction  between  a  shock  and  body  geometry  should 
not  be  essential,  so  long  as  the  shock-layer  approximation  is  valid,  since 
the  two  are  equivalent  in  the  leading  approximation.  To  provide  a  flexible 
interpretation  of  the  present  analysis,  one  may  regard  either  the  body 
shape  or  the  shape  of  the  "shock  interface"  as  specified. 


The  following  study  concerns  highly  cooled  bodies  of  uniform  surface 

temperature  with  ,  8  —  and  Pf  -711  ■  To  simplify  the 

I  o  o  0 

matter,  the  linear  viscosity-temperature  relation  —  ~r/7^  is 

assumed.  The  constant  / 7^  can  be  evaluated  in  the  manner  des- 


ican  be  evaluated  in  the  manner  des¬ 


cribed  in  Section  4.  1  of  Part  I,  using  the  reference  temperature  given  in 


Eq.  (4.  la)  for  the  stagnation  region. 


This  linear  representation  may  perhaps  be  somewhat  improved  by  taking 
into  account  the  effect  of  flow  speed  on  the  mean  temperature  of  the  shock 
layer,  in  a  manner  similar  to  that  provided  in  Eq.  (7.4). 
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In  the  present  study  concerning  flow  regions  not  very  far  from  the 
blunt  nose,  such  idealization  should  be  sufficient.  In  terms  of  the  variable  A  , 
the  shock  or  body  shape  considered  can  be  described  as 


(7.8) 


where  fo  is  zero  for  a  paraboloid,  and  is  positive  for  a  hyperboloid.  The 
analyses  to  be  discussed  are  obtained  for  b  =  0  and  b  =  1  .  The  latter 
corresponds  to  a  hyperboloid  which  approaches  a  45°  cone  at  infinity  -- 
(ideally  the  hyperboloid  may  be  regarded  as  a  blunted  cone).  For  each 
geometry,  values  of  )<C  -  0.3,  1,  3,  10  and  30  are  considered,  encom¬ 

passing  the  viscous-layer,  the  incipient-merged  layer  as  well  as  the  fully- 
merged  layer  regimes. 

The  method  of  series  expansion  is  first  employed  to  determine  Cce  ,  ©0  , 

-f2',  >  Uf  <  0/  >  an8  which  provides  initial  data  for  the  difference 

2 

method.  The  solutions  are  carried  out  for  each  value  of  «  from  A  =0.1 
near  the  axis  of  symmetry  to  A  =  3  (about  2  to  3  times  the  nose  radius), 
with  a  step  size  of  AX  =  1/100  and  =  1/100.  To  avoid  exce  s  sive 

amplification  of  the  truncation  error  near  ^  =  0,  the  inner  boundary  is 

shifted  from  £  =  0  to  £  =  ,  as  previously  discussed.  The  value  of  C, g 

used  is  0.15  for  the  cases  then  l<  =0.3,  1,  and  3,  and  0.10  for  the  less 
critical  cases  where  ^  =10  and  30.  Instability  is  not  detected,  except 

for  a  bounded  oscillation  of  small  amplitude  which  is  found  in  the  case  of 
K  =  0.  30  over  the  range  of  0.  1  £  A&  0.  30.  For  each  value  of  K*  ,  the 
total  computation  time  (including  print  out)  on  an  IBM  704  digital  computer 
is  about  27  minutes.  The  results  of  the  calculation  are  presented  in  Figs. 

7.3  to  7.11. 
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The  profiles  of  velocity,  enthalpy,  pressure  and  temperature  at  three 
successive  stations  —  =  .59,  1.24  and  2.82  downstream  of  the  stagnation 
point  are  presented  in  Figs.  7.3  to  7.5  for  various  values  of  K  The 
profiles  for  the  paraboloidal  case  which  are  essentially  similar  are  not 
presented.  The  particular  features  of  these  profiles  which  should  be 
emphasized  are  the  reductions  of  tangential  velocity  and  total  enthalpy  at 
the  outer  edge  from  their  respective  upstream  values  as  a  result  of  the 
transport  effects  behind  the  shock.  The  reductions  increase,  of  course, 
with  decreasing  density,  i.e.,  decreasing  «■  ,  and  decrease  generally 

with  increasing  distance  rC  from  the  axis  of  symmetry.  As  in  the  solution 

,  etc . 

near  the  body  surface  increase  with  distance  ^  ,  and  a  boundary  layer  is 

seen  to  emerge  gradually  at  the  base  of  the  shock  layer,  as  one  moves 

downstream.  At  the  higher  Reynolds  numbers,  K"  *  =  1  0  and  30,  the 

temperature  profiles  at  the  downstream  stations  attain  maximum  values 

inside  the  shock  layer,  as  in  a  boundary  layer  with  surface  cooling.  The 

profile  of  the  total -enthalpy  function  ©  also  reveals  similar  characteristics, 

though  to  a  lesser  degree.  In  fact,  for  the  case  of  l(  =  30,  the  function  &  , 

thus  the  total  enthalpy  itself,  overshoots  slightly  its  free- stream  value  at 

the  station  -  =  2.82.  This  overshoot  results  from  the  mutual  compensa¬ 

te 

tion  of  the  viscous  dissipation  and  heat  conduction  (from  the  higher- 
temperature  region  inside  the  shock  layer),  with  the  former  overbalanced 
by  the  latter.  This  feature  is  peculiar  only  for  the  Prandtl  number  0.  71 
assumed,  and  is  absent  for  Pr  >  |  .  The  similar  overshoot  in  total  enthalpy 
has  also  been  found  previously  in  Van  Driest's  calculations  of  laminar  boundary 


for  shock  layers  on  the  wedge  and  cone, 


'lauicma  OU. 


'  J 
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layer  over  a  flat  plate.  One  may  also  note  from  these  results  that  there 
is  very  little  variation  of  pressure  across  the  shock  layer  for  all  Reynolds 
numbers.  This  property  is,  of  course,  only  incidental  to  the  particular 
configuration  studied. 

From  the  velocity,  temperature  and  pressure,  the  streamline  pattern 
describing  the  flow  field  around  a  blunt  body  can  be  determined.  Results  are 
presented  in  Figs.  7. 6  and  7. 7  for  «  =1  and  30  for  both  the  paraboloid 

and  the  hyperboloid.  In  order  to  provide  a  closer  comparison  with  the 
inviscid  case  previously  presented  in  Fig.  5.2  of  Part  I,  the  results 
presented  in  Figs.  7. 6  and  7.  7  have  been  interpreted  as  solutions  to  the 
inverse  problem,  that  is,  the  flow  fields  and  body  surface  supporting 
paraboloidal  and  hyperboloidal  shocks.  To  be  more  specific,  by  the  term 
"shock"  one  is  referring  to  the  "shock  interface"  discussed  in  Section  5, 
which  is  also  the  location  of  the  outer  edge  of  the  shock  layer  in  the  present 
two-layer  model.  Consider  first  Fig.  7 . 6  at  K  =1.0,  which  represents 
the  flow  fields  in  the  merged-layer  regimes.  Comparing  the  results  of  the 
paraboloidal  shock  in  Fig.  7.6(a)  with  the  corresponding  results  in  Fig.  5.2, 
one  sees  that  far  from  the  axis  of  symmetry  the  inviscid  and  viscous  shock 
layers  are  comparable  in  thickness,  but  that,  near  the  stagnation  region, 
the  thickness  in  the  viscous  case  studied  is  noticeably  smaller.  The  results 
are  understandable  through  the  observation  that  the  high  rate  of  surface 
cooling  effectively  increases  the  density  level,  and  hence  the  mass  flux, 
p  ,  near  the  stagnation  region.  Away  from  the  stagnation  region,  the 

Strictly  speaking,  the  particular  example  given  in  Fig.  5.2  which  assumes 
£  =1/7  does  not  correspond  to  the  inviscid  limit  of  the  example  considered 

here,  which  assumes  £  =  1/8.  The  difference  however  is  rather  small. 
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heat  conduction  to  the  cooled  body  is  offset  and  even  overbalanced  by  the 
viscous  dissipation.  For  an  insulated  body,  one  may  anticipate  a  sub¬ 
stantial  thickening  of  the  shock  layer  downstream  of  the  stagnation  point. 
Comparing  the  flow  pattern  of  the  paraboloidal  (Fig.  7.  6a),  with  the  hyper  - 
boloidal  shocks  (Fig.  7.6b),  the  latter  is  seen  to  be  generally  thinner.  The 
relative  difference  between  the  inviscid  and  the  "fully  viscous"  shock  layers 

observed  for  the  paraboloid  case  also  applies  for  the  hyperboloidal  case. 

2. 

The  flow  patterns  at  a  higher  Reynolds  number,  K  =  30,  is  pre¬ 
sented  in  Fig.  7.  7.  It  is  of  interest  to  observe  that  the  shock-layer  thicknesses 
around  the  bodies  have  not  changed  much  from  those  of  the  preceding  figure, 
in  spite  of  the  considerably  larger  difference  in  Reynolds  number.  The  main 
reason  is  presumably  that,  at  -  30,  a  good  fraction  of  the  shock  layer 
near  the  stagnation  point  still  belongs  to  the  region  where  viscous  and  heat- 
conducting  effects  are  important. 

It  may  be  remarked  that  the  re  suits  pre  sented  in  Figs .  7.6  and  7 . 7, 
as  well  as  in  the  other  figures,  do  not  provide  the  complete  flow  picture, 
since  they  describe  only  the  pattern  within  the  shock  layer.  As  pointed 
out  before,  solution  to  the  shock-transition  zone  can  be  easily  obtained  from 
the  information  provided  here.  The  streamline  pattern  in  the  shock- 
transition  zone  may  be  constructed  by  working  backward. 

The  surface  pressures  on  the  paraboloid  and  hyperboloid,  non- 
dimensionalized  by  the  pressure  behind  the  shock,  are  provided  in  Fig,  7.8. 
Although  not  corresponding  to  the  exact  inviscid  limit  for  £  =  1  /8,  the 
inviscid  results  calculated  previously  for  £  =1/7  are  also  included  for 
comparison.  The  results  obtained  show  a  continuous  transition  between  the 
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inviscid  limit  and  fr ee-molecule  limit  (which  is  ■P  =  1}-  AS  already 
noted,  i2'  is  practically  unity  for  the  hyperboloid.  Also  included  in  Fig. 

7. 8  are  the  distribution  of  the  thickness  of  the  shock  layer  around  bodies 

2,  2. 
for  the  five  cases  of  K  analyzed.  Except  for  the  case  of  K  =  0.  30,  the 

results  for  all  K  *  1  s  compare  closely  one  with  another,  confirming  the 

relative  insensitivity  of  the  shock-layer  thickness  with  respect  to  change  in 

Reynolds  number  previously  observed. 

Finally,  the  local  skin  friction  and  surface  heat-transfer  rate  are 

presented  respectively  in  Figs.  7.8  and  7.9  for  both  the  paraboloid  and 

hyperboloid.  Also  included  are  the  distributions  corresponding  to  the 

free-molecule  limit  (for  unit  accommodation  coefficients).  The  smooth 

transition  to  the  free-molecule  limit  is  quite  evident  from  the  results 

x 

presented  for  successive  values  of  fC  .  Generally,  the  value  of  <ZH  falls 
off  as  one  moves  away  from  the  stagnation  point.  Similar  is  the  value  of  , 
after  reaching  a  maximum.  One  notes  that  the  rate  of  this  falling-off  increase 
generally  with  increasing  Reynolds  number,  also  that  the  distributions  of 
the  higher  's  arc  remarkably  similar.  To  study  the  heat-transfer  charac¬ 
teristics  more  closely  and  to  make  comparison  with  the  boundary-layer 
limit  possible,  the  results  of  Fig.  7.  10  are  reproduced  in  Fig.  7.11  in  the 
form  of  ^■(V,Ax  vs.  /<*■  ,  where  is  the  heat-transfer  rate  at 

the  stagnation  point.  One  sees  that  for  both  paraboloid  and  hyperboloid,  the 
heat-transfer  distributions  at  kC*"3  30,  10,  as  well  as  3,  are  reasonably 
close  to  one  another.  While  the  rate  of  falling-off  from  the  stagnation-point 
value  is  seen  to  increase  with  increasing  K  as  in  the  preceding  figure, 

Fig.  7.11  reveals  that  there  is  a  reversal  of  this  trend  for  at 
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«  =10  for  both  paraboloid  and  hyperboloid.  That  is,  at  the  same  qt/a.  , 

the  ratio  a. / is  slightly  higher  for  K%=  30  than  for  kT  *  =  10.  This 
reversed  trend  is  further  confirmed  by  the  boundary-layer  limit  estimated 
on  the  basis  of  the  local- similarity  method  (using  the  inviscid  solution  ob¬ 
tained  from  the  present  theory  for  the  outer  flow),  also  included  in  Fig.  7.11. 

2. 

The  reversal  in  the  trend  of  heat-transfer  distribution  at  K  =10  can  be  ex¬ 
plained  by  the  vorticity  effect.  Since  the  external  vorticity  effect  generally 
decreases  in  the  downstream  direction,  the  effect  may  tend  to  decrease  on 

ratio  from  the  corresponding  boundary- layer  prediction,  in  agreement 
with  the  reversed  trend  observed. 

From  the  above  discussion,  one  should  note  that  the  distribution  of  the 

heat-tr ansfer  ratio  is  relatively  independent  of  Reynolds  number  over  a  very 

wide  range.  In  particular,  the  variation  of  the  heat-transf er  ratio  in  the  range 
2. 

K  =  3-30  is  remarkably  small.  This,  together  with  the  fact  that  the  vorticity 
interaction  is  itself  not  a  predominantly  large  effect,  suggests  that  the  heat- 
transfer  distributions  on  a  smooth  blunt  body  should  remain  essentially  un¬ 
changed  in  the  boundary-layer  through  the  incipient-merged-layer  regimes. 

A  similar  conclusion  is  applicable  to  the  skin-friction  characteristics. 

The  insensitivity  of  the  distribution  of  the  measured  heat-transfer  rate  with 
respect  to  Reynolds  number  has  been  reported  in  Refs.  50  and  52  for  low- 
density  hypersonic  flows  over  spheres. 


The  accuracy  of  the  local-similarity  method  for  the  boundary-layer  analysis 
is  taken  for  granted  here. 


In  this  instance,  one  shall  consider  the  distribution  of  the  ratio  of  f  to 

J&m,  /%,  d.f . 
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8.  STUDY  OF  NONEQUILIBRIUM  DISSOCIATION  IN  SHOCK  LAYER  AT  LOW 
REYNOLDS  NUMBER 

Thus  far,  only  an  ideal  gas  with  constant  specific  heats  has  been  considered. 
A  more  realistic  analysis  of  the  low-density  hypersonic  flow  must  include  non¬ 
equilibrium  flow  chemistry.  The  following  section  will  discuss  briefly  an  exten¬ 
sion  of  the  thin-layer  approach  to  study  dissociation  and  species  diffusion  in  the 
viscous  shock  layer,  assuming  a  single  dissociation-recombination  reaction. 

The  stagnation  region  in  the  viscous  layer  and  merged-layer  regimes  will  be 
analyzed. 

8.1  Relation  to  Previous  Work 

The  stagnation  region  with  nonequilibrium  dissociation  has  been  studied 
3  6 

by  Chung,  on  the  basis  of  a  viscous-layer  model.  His  analysis  assumes 
a  binary  mixture  of  diatomic  molecules  with  a  unit  Lewis  number  and  a  linear 
viscosity-temperature  relation,  for  a  noncatalytic  surface.  Very  significant 
effects  of  surface  cooling  on  the  flow  chemistry  have  been  revealed.  The 
following  analysis  assumes  a  similar  chemical  kinetic  model  but  will  be  con¬ 
ducted  from  a  viewpoint  more  consistent  with  the  general  approach  developed 
in  the  preceding  sections.  The  present  flow  model  allows  for  changes  in  the 
tangential  velocity,  total-enthalpy,  and  species  concentration  across  the  shock 
and  therefore  possesses  a  regime  of  validity  larger  than  the  viscous -layer 
model.  The  influence  of  the  wall  catalycity  on  the  species  diffusion  and  flow 
chemistry  will  also  be  studied. 

In  considering  nonequilibrium  flow  of  rarefied  gases,  questions  naturally 

arise  as  to  the  importance  of  the  effects  associated  with  the  vibrational  (also 

79 

perhaps  rotational)  relaxation.  Although  not  being  accounted  for  in  the 
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present  study,  the  basic  flow  model  developed  here  may  provide  a  simple 
framework  for  studying  these  as  well  as  other  more  complex  high-temperature 
real-gas  effects. 

Before  passing  to  the  specific  analysis  of  the  stagnation  region,  the 
general  problem  of  the  shock  layer  with  species  diffusion  and  nonequilibrium 
chemistry  will  be  discussed. 

8.2  General  Discussion 


Consider  a  binary  mixture  of  diatomic  molecules  and  dissociated  atoms. 

The  diffusion  velocity  of  the  atomic  species  V '  may  be  written  for  the  present 

,  ,  6,80 
.  study  as 

v/=-D'*  "(fr^+oC)  (8.1 


where  oc  is  the  atomic  concentration  in  mass  fraction,  and  £>  l2  the  binary 
diffusion  coefficient.  The  thermal  diffusion  has  been  generally  neglected  in 
boundary  layer  analysis  (  1  / 1 0  for  air)  but  may  be  quite  significant 

for  nonequilibrium  flow  with  a  highly  cooled  surface. 

Of  particular  significance  is  the  second  term  on  the  right  of  Eq.  (8.1) 
corresponding  to  the  pressure  diffusion.  This  term  is  not  present  in  boundary- 
layer  analysis  which  assumes  a  uniform  pressure  across  the  layer,  but  is 
important  in  the  present  study  because  the  pressure  variation  is  not  generally 
negligible  in  the  shock  layer  (as  well  as  in  the  shock-transition  zone).  In  fact, 
in  the  merged-layer  regimes,  the  contribution  of  the  pressure  diffusion  can 
be  as  important  as  that  due  to  any  other  processes,  and  should  therefore  be 
included  in  the  general  formulation  of  the  problem.  In  the  stagnation  region, 
however,  the  pressure  gradient  is  small  and  the  pressure  diffusion  effect  can 
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also  be  neglected.  From  Eq.  (8.1),  one  may  infer  that  pressure  diffusion  tends 
to  increase  the  flux  of  the  (lighter)  atomic  species  toward  the  (convex)  surface 


in  region  downstream  of  the  stagnation  point. 

6  V  V 

Consistent  with  Eq.  (8.1)  is  the  expression  for  the  energy  flux  ’ 

t  =  "B  IS  u  (8.2) 

6  a 

where  Q  is  the  heat  of  formation  per  unit  mass  of  the  atoms,  and  and 
the  specific  enthalpies  of  the  atoms  and  molecules,  respectively.  The  differ¬ 
ential  equations  governing  the  shock  layer  Eq.  (5.3)  remain  unchanged  except 
for  the  last  equation  governing  the  total  enthalpy.  The  equations  governing  the 
total  enthalpy  H  and  the  atomic  concentration  oC  in  the  shock-layer  region 
are  i-  — . 

p(u  =  sf- 

where  M)~'  is  the  rate  of  production  of  the  atomic  species  in  mass  per  unit 
volume.  Consistent  wir.h  the  present  two-layer  model,  and  with  Eqs.  (8.1) 
and  (8.  2),  are  the  two  equations  relating  total -enthalpy  drop  and  concentration 

/c 

The  assumption  of  and  Q  has  been  used. 

To  complete  the  system  of  differential  equation,  one  has  -p- R2(l  +  cc)pT , 
where  is  the  gas  constant  for  the  molecule,  as  well  as  X,  and  as 
function  of  T  ,  /(  as  function  of  T  and  o i  ,  etc. 
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jump  across  the  shock-transition  zone. 

p,  v,  ('oc2-t>;00)=  -(<*  v')z 

The  finite  atom  concentration  at  the  outer  edge  of  the  shock  layer,  ,  is 
a  result  of  upstream  diffusion  of  the  atomic  species  which  are  produced  inside 
the  shock  layer.  Implicit  in  the  last  of  Eqs.  (8.4),  which  resemble  a  concentration- 
jump  across  the  Knudsen  layer,  is  the  assumption  that  the  gas-phase  reaction  is 
frozen  within  the  shock-transition  zone,  i.e.  aas*  is  negligible.  The  assumption 
may  be  justified  on  the  basis  of  the  relatively  low  collision  frequency  in  the 
transition  zone.  Because  of  the  temperature  maximum  reached  inside  the  trans¬ 
ition  zone  however,  a  domain  (corresponding  to  a  flight  speed  and  altitude  much 
higher  than  contemplated)  does  exist  in  which  excitations  of  modes  of  vibration 
and  dissociation  within  the  shock  can  be  important. 

The  boundary  condition  governing  the  species  concentration  at  the  surface, 

78 

assuming  a  first-order  recombination  rate  and  neglecting  dissociation  at  the 
surface,  is 


~(ocV')  = 


oC 


y 


(y  -  o') 


(8.5) 


where  '/  is  the  recombination  efficiency  (the  probability  for  each  atom  to 

w 

recombine  after  reading  the  surface)  and  ^0C\/'_)is  provided  by  Eq .  (8.1)  in¬ 
cluding  the  pressure  diffusion.  To  see  more  clearly  the  manner  in  which  the 
recombination  efficiency  affects  surface  catalycity  in  the  regimes  con¬ 
sidered,  one  may  express  Eq.  (8.5)  in  terms  of  the  variable  £  , 

/+y  /z-rw\.r  T7  A,,  *  a* 


oC  = 


ZK;  Pr 

o 


U rj7  ^  £ 


0$ 


(8.6) 
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One  sees  that  as  Kfl  decreases,  for  a  fixed  recombination  efficiency 


the  wall-catalycity  effect  reduces.  In  the  range  of  =  0(1),  the  question  of 

whether  or  not  a  fully  catalytic  wall  (i.e.  OC  O  )  can  be  realized  will  depend 
on  whether  _ 

I  LL  ,  /t~ 

«  / 


I  _A 
Ac 


w 


/ 


R 

T 


w 


(8.  6a) 


Since  /<  °C  ,  and  cJ  is  not  much  larger  than  1/2,  it  seems  that  a  fully 

catalytic  surface  would  require  comparatively  high  value  of  7^  at  =  0(1). 

Assuming  a  single  dissociation- recombination  reaction  of  the  type 
+  X  +  X+M,  the  rate  of  production  Atf/p  may  be  written  in  the  familiar 


form 


28, 29,79 

~ 


(8.7) 


where  is  the  "equilibrium  constant"  which  may  be  related  to  the  constant 

WAT 


19 

p>  of  Lighthill  through 


(8.8) 


P  =  J£p_f_ 

+RtT 

Lighthill  finds  p  remains  practically  constant  for  pure  oxygen  and  pure 

nitrogen  for  a  temperature  as  high  as  8000°K.  Examination  shows  that  at 

higher  temperature,  large  departure  of  p  from  the  nearly  constant  value 

will  occur  for  nitrogen  because  of  the  excitation  of  certain  highly  degenerated 

electronic  states  of  the  nitrogen  atoms.  In  his  study  of  the  viscous  shock  layer, 

36 

Chung  assumed  "unexcited  electronic  state"  (apparently  meaning  the  high- 
lying  levels  to  be  unexcited).  Presumably,  his  p^  does  not  differ  a  great 
deal  from  its  value  at  8000°K. 

Another  point  of  importance  with  regard  to  the  reaction  rate  is  the  re¬ 
latively  weak  recombination  in  the  gas  phase  at  low  density.  As  observed  by 
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84 

Gibson,  the  recombination  rate  based  on  the  three-body  process  becomes 

29  30  82 

negligible  as  density  level  decreases.  Inboundary-layer  studies,  ’  ’ 

the  recombination  is  important  because  the  ambient  density  level  is  sufficiently 
high  and  increases  further  towards  a  highly  cooled  surface.  In  the  low  Rey- 
olds  number  regimes  studied  here,  the  ambient  density  is  too  low  to  make  the 
three-body  processes  effective,  even  with  a  highly  cooled  surface.  The  re¬ 
combination  in  the  gas  phase  may  therefore  be  expected  to  be  negligible.  The 
range  of  validity  of  this  simplification  can  be  easily  established  by  an  order-of- 
magnitude  analysis  based  on  Eq.  (8.  7).  A  direct  consequence  of  neglecting  the 
three-body  processes  is  the  binary  scaling  law,  which  permits  correlations 
of  nonequilibrium  flow  fields  of  different  body  scale  and  ambient  density  by 

preserving  the  product  (  pa,  )  in  the  same  manner  as  in  the  Reynolds-number 
84 

similitude.  It  is  evident  that  this  scaling  law  may  apply  to  both  inviscid  and 
viscous  flows.  In  case  of  viscous  flow,  however,  in  order  to  preserve  similar¬ 
ity  in  the  boundary  condition,  the  binary  scaling  law  requires  the  nonslip  surface 

doc 


condition,  and  a  noncatalytic  (  — •  =  0)  or  a  fully  catalytic  surface  (  oC  =  0). 

In  the  following,  nonequilibrium  dissociation  in  the  stagnation  region  will  be 
analyzed  on  the  basis  of  the  model  described.  In  passing,  it  may  be  pointed 
out  that  the  procedure  of  solution  developed  in  the  previous  section  for  shock 
layers  around  nonslender  bodies  can  be  extended  to  study  nonequilibrium 
dissociation,  including  species  diffusion  in  the  shock-transition  zone. 

8.3  Nonequilibrium  Dissociation  and  Species  Diffusion  in  the  Stagnation  Region 
In  the  following  study,  the  perfect-gas  value  of  for  the  atom  species 
and  molecular  species  will  be  assumed  to  be,  respectively, 
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=i  *. 

cr.=(r+')«.=fe*  J  (8 

Thus  the  molecules  in  the  shock  layer  are  assumed  to  possess,  in  a  classical 
sense,  the  full  rotational  as  well  as  vibrational  degrees  of  freedom.  Implicit 
is  the  assumption  that  excitation  of  the  electronic  states  of  the  atoms  and 
molecules  does  not  contribute  to  the  internal  energy  of  the  particle.  Thus, 
associated  with  the  atomic  and  molecular  species  are  two  specific  heat  ratios 
■y  =  5/3  and  y  -  9/7.  To  make  the  temperatures  in  the  following  analysis 


dimensionless,  one  introduces  a  "free-stream  stagnation  temperature": 


T  = 


(8.1i; 


where  corresponds  to  a  specific  ratio  of  7^  =7/5. 

00 

It  will  be  convenient  to  introduce  the  following  parameters 

e*  3  V-iT.  'j 

*  =  2 V 
c<  u‘r“ 

0O  4  00 

<r  =  p  A  p/  J 

'00/  00  I  d 

where  is  a  reference  value  for  p  and  will  be  taken  for  the  preset 

purpose  as  122  gm/cc. 


(8.12) 


The  system  of  differential  equations  governing  the  shock  layer  in  the 


stagnation  region  in  terms  of  the  variable  C,  ,  becomes 
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where  H  —  H/^  and  &^T/T0  are  respectively  the  dimension^ 
enthalpy  and  temperature,  which  are  related  through 

97t(-T*-e*«) 

iat  K  (*6  *■)  is 


(8.14) 


Note  that 


is  the  same  parameter  as  defined  previously.  In 


the  above  equations,  a  viscosity  law  JU.  °C  T^has  been  used,  and  the  thermal 
diffusion  neglected. 

The  boundary  condition  on  the  surface  becomes,  at  £  =  0 

a  =0 

e  =  VTo 


(8.15) 


where 


06  "  U.'oc'  =0 


r  = 


/  v4*-* 

Zrw  /  05  VTo/ 


The  boundary  condition  at  the  outer  edge  of  the  shock  layer  becomes,  at  =  1 


a  + 


u_u' 


C  e'"w(/  +  °0 


ft  =0, 


..  .  <*fe  Hoc' 


-  Oc 


(8.16) 


H  + 


Pr*C 


ft  H  S«o(cJfe  -iXS  +  ^9^)  ft  0^  _  1  L  o  ^  /n 
Pr  IT,1  ©'-»(*  +  «) - ' 


*  co  ^  ✓ 

Eqs.  (8.13)  to  (8.  16)  constitute  a  standard  two-point  boundary- value  problem. 
The  above  formulation  includes  the  frozen  atomic  species  in  the  free  stream. 

In  the  following  examples,  the  above  formulation  is  applied  to  an  ideal¬ 
ized  mixture  of  "air  molecules"  and  air  "atom"  as  in  Refs.  29  and  79,  thus 
one  uses  a  value  of  =  T>/&  =  1  01 , 300°  K,  which  is  a  weighted  average 

between  59,000°K  (oxygen)  and  113,  200°  K  (nitrogen),  based  on  the  composition 


This  idealization  may  lead  to  oversimplification  of  the  chemical  rate  equation. 
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in  the  free  stream.  To  make  the  present  calculations  comparable  with  Chung's 
analysis,  the  same  values  of  (J&  ,  ZJk'  ,  S  ,  a.  ,  and  oC^  (=0),  are 

used  and  are  listed  in  Fig.  8.  1.  Differing  from  Chung's  however  are  the 
values  CO  =1/2  and  dS=  =  1.40,  also  p  Ip1  -  1.  These  differences  should 
not  be  considered  to  be  critical.  The  value  U ^  =  26,  000  ft/sec  used  corres¬ 
ponds  to  re-entry  from  satellite  speed. 

Numerical  integrations  have  been  performed  for  =  0.  473,  1 . 97, 

5. 79,  16.  24  and  35.2,  corresponding  to  a  one-foot  nose  radius  at  altitudes 
of  350,  000  ft.,  320,000  ft.,  300,000  ft.,  280,000  ft.  and  265,000  ft., 


respectively.  Throughout  the  whole  series  of  calculations,  the  values  of 
\<y^QZ+^  which  controls  the  recombination  rate  are  of  the  order  of  1  0  2  or 
much  less,  and  one  may  use  the  binary- scaling  rule  to  apply  the  present 

results  to  bodies  of  other  dimension.  The  table  given  on  top  of  Fig.  8.  1 

z 

provides  the  correspondence  between  the  value  of  KQ  and  the  altitude  for 

nose  radius  smaller  and  greater  than  one  foot,  according  to  this  rule. 

The  solutions  by  numerical  integration  are  presented  in  Fig.  8.  1  to 

8.3.  As  typical  results,  the  atom  concentration  profiles  of  the  shock  layer 

are  presented  in  Fig.  8.  1  for  both  noncatalytic  as  well  as  fully  catalytic 
2  2. 

surfaces  for  =  16.24.  This  value  of  kCo  falls  in  the  viscous-layer 

regime  and  corresponds  to  the  highest  altitude  considered  by  Chung  (280,000 

ft.).  Included  in  Fig.  8.  1  for  comparison  is  the  corresponding  profile  in  the 

2 

inviscid  limit  (  K0 — >oo),  as  well  as  a  solution  for  unit  Lewis  number.  The 
difference  in  solution  due  to  the  Lewis  number  is  clearly  small,  at  least 
for  the  noncatalytic  case  considered,  whereas  the  departure  of  the  solutions 
from  the  inviscid  limit  is  clearly  very  large. 
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An  important  feature  of  the  profiles  obtained  is  the  finite  value  of  atom 
concentration  at  the  outer  edge,  which  is  absent  in  the  corresponding  solution 
of  Chung.  The  present  results  show  that  the  concentration  jump  is  not  neg¬ 
ligible  even  in  the  viscous-layer  regime  where  changes  in  tangential  velocity 
and  enthalpy  across  the  shock  may  be  neglected.  Chung’s  solution  (not  shown) 
follows  the  inviscid  profile  near  the  shock  and  becomes  comparable  with  the 
present  result  (for  the  noncatalytic  case)  near  the  surface.  The  large  influence 
of  the  wall  catalycity  on  the  composition  in  the  inner  portion  of  the  shock  layer 
is  clearly  indicated  by  the  large  difference  between  the  two  extreme  cases 
shown  in  Fig.  8.1.  However,  the  temperature  profiles  corresponding  to  these 
extreme  cases  (not  shown)  reveal  indistinguishable  difference  throughout  the 
entire  shock  layer.  The  negligible  difference  in  temperature  profile  is  ex- 
plained  by  the  fact  that  for  this  as  well  as  other  values  of  considered 

(and  for  the  speed  of  23,  000  ft/sec  or  lower)  dissociation  takes  place  mainly 
in  the  outer  portion  of  the  shock  layer  where  the  temperature  is  highest.  The 
composition  in  the  inner  portion  of  the  shock  layer  is  controlled  by  diffusion 
as  in  a  frozen  boundary  layer.  Since  there  is  practically  no  heat  release  from 
chemical  reaction  in  the  inner  part  of  the  layer,  the  change  in  composition 
resulting  from  the  change  inwall  catalycity  will  not  introduce  appreciable 
change  in  temperature. 

The  profiles  of  concentration,  velocity,  enthalpy  and  temperature  at 
2 

other  values  of  (<0  corresponding  to  successive  degrees  of  rarefaction  are 
presented  in  Fig.  8.  2  for  the  case  of  noncatalytic  surface.  The  corresponding 
results  for  a  fully  catalytic  surface  will  not  be  presented,  except  to  mention 
that  the  conclusion  drawn  from  the  study  of  Fig.  8.  1  with  respect  to  the  large 
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influence  of  the  wall  catalycity  on  atom  concentration  profile  and  the  small 
effect  on  temperature  remains  the  same.  According  to  Fig.  8.2,  at  a 

%  i 

higher  value  of  K0  ,  say  Ko  =  35.2,  the  concentration  jump  across  the 
shock  is  still  appreciable  and  reduction  of  atom  concentration  level  by 

a  2 

surface  cooling  is  still  significant.  At  the  lower  values  of  ,  say  Ka  —  5, 

the  temperature  is  reduced  (because  of  the  combined  effect  of  surface  cooling 
and  enthalpy  drop  across  the  shock)  to  such  a  low  level  that  the  effect  of  dis¬ 
sociation  is  hardly  distinguishable. 

The  stagnation-point  heat-transfer  coefficient  C  H  as  well  as  the 
standoff  distance,  A  ,  corresponding  to  the  solutions  discussed  are  presented 
in  Fig.  8.  3.  As  anticipated,  the  use  of  a  non-catalytic  surface  to  reduce  heat 

transfer  becomes  less  and  less  effective  as  the  degree  of  rarefaction  increases, 

z 

i,  e.  K0  decreases,  since  the  difference  in  heat-transfer  characteristics  of 

the  two  extreme  surface  conditions  depend  on  the  amount  of  atomic  particles 
2 

present.  At  Ka  below  5,  the  difference  in  heat  transfer  vanishes,  Also  in¬ 
cluded  is  the  CH  value  for  a  Y  =  9/7,  corresponding  to  the  case  of  completely 
frozen  chemistry.  The  stand-off  distance,  that  is,  the  thickness  of  the  shock 
layer  (obtained  for  a  noncatalytic  surface)  is  seen  to  vary  with  density  rather 
slowly.  The  relatively  high  degree  of  dissociation  taking  place  in  the  outer 
portion  of  the  shock  layer  at  higher  «0  ‘s  reduces  the  temperature  and  in¬ 
creases  the  density,  causing  a  reduction  in  the  standoff  distance  from  the 
ideal-gas  value,  as  shown  in  Fig.  8.  3. 
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CONCLUSION 


In  this  paper,  the  problem  of  hypersonic  flow  of  a  rarefied  gas  around 
blunt  bodies  is  studied  within  the  framework  of  the  continuum  theory.  The 
existing  theories  and  experiments  on  the  stagnation  region  are  discussed;  the 
approach  based  on  the  thin- shock-layer  approximation  is  extended  to  study 
flow  fields  removed  from  the  stagnation  region  as  well  as  in  the  shock- 
transition  zones.  Particular  problems  are  analyzed  to  illustrate  aerodynamic 
and  heat-transfer  characteristics  of  blunt  as  well  as  pointed  nonslender  bodies 
at  low  Reynolds  number.  Also  studied  is  nonequilibrium  di ssociation  with 
species  diffusion  in  the  shock  layer  under  the  influence  of  strong  surface 
cooling. 

Attempts  to  analyze  the  stagnation  region  based  on  the  continuum  theory 

have  resulted  in  three  main  approaches.  One  approach  treats  vorticity  and 

other  effects  under  the  assumption  of  small  departure  from  the  boundary- 

24  25  26 

layer  limit.  The  theories  of  Van  Dyke  and  others  ’  based  on  this  approach 
have  been  most  systematically  and  elegantly  developed,  and,  in  fact,  can  be 
extended  to  study  flow  regions  removed  from  the  stagnation  point.  In  a  second 
approach,  applicable  to  regimes  of  low  density,  transport  effects  are  considered 
throughout  the  shock  layer,  and  a  system  of  simplified  Navier-Stokes  equations 
is  integrated  along  the  axis  of  symmetry.  However,  in  order  to  make  the 
partial  differential  equations  tractable  in  this  approach,  one  has  to  assume  a 
constant  density  or  a  thin  shock  layer.  The  third  approach  makes  consistent 
use  of  the  thin- shock-layer  model.  Without  sacrificing  the  essential  features 
of  the  problem,  the  thin- shock-layer  approach  simplified  the  analyses  of  the 
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stagnation  region  in  both  the  high  and  low  Reynolds  number  regimes 
Within  the  degree  of  approximation  used,  it  embraces  both  of  the  other 
approaches 

As  a  result  ot  the  more  systematic  studies,  ^  ^  much  ol'  the  contro¬ 
versy  and  disagreement  in  results  of  vor ticity ,  displacement  and  other  eifects 
has  bt-n  resolved  during  recent  years,  unresolved,  however,  remains  the 
variance  between  the  predicted  and  measured  values  of  heat-transfer  rates 

One  may  cite  in  pa rticular  that  the  rise  01  heat-transfer  rate  predicted  by 
26  24  45 

Maslen  Van  Dyke  and  others  is  consistently  cower  than  the  values 


deduced  from  the  experimental  measurement  of  Ferri,  Zakkay  and  Ting, 


27 


and  some  other  sources  Unresolved  also  remains  the  difference  among  results 

2  6 

tor  the  slip  and  temperature-jump  effects,  Maslen  shows  that  these  eifects 

are  net  important  for  heat-transfer  studies  as  opposed  to  the  conclusion  oi 

Rotr  and  Len^rd  ^  and  also  of  Van  Dyke.  ^  Disagreement  appears  to  be 

Its-  critical  among  results  for  the  lower  Reynolds  number  regime.  Experi- 

*4  7 

mental  heat -transfer  data  measured  by  Ferri  and  Zakkay.  and  by  Vidal  and 

Wittliff^  appear  to  follow  the  theoretical  trend  of  the  thin- shock-layer  analysis 

rather  weii  in  this  regime.  There  exists  also  reasonable  correlation  among 

heat-transfer  predictions  of  Probstein  and  Kemp,  J  ^  Levinsky  and  Yoshihara"^ 

41 

and  Cheng  in  the  merged-layer  regimes, 

The  new  development  with  the  thin- shock-layer  approach  includes  an 
extension  of  the  thin-layer  concept  to  formulate  the  problem  in  the  shock  - 
transition  zone.,  and  development  of  a  finite -difference  method  for  solving  the 
shock-layer  equations  without  similarity  assumptions 

The  basic  flow  model  consists  of  an  outer  layer  and  an  inner  layer 
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identifiable  as  the  shock-transition  zone  and  the  shock  layer,  respectively. 
Across  the  shock-transition  zone,  the  compression  ratio  is  assumed  to  be 
high.  In  the  shock- transition  zone,  the  density  is  assumed  to  be  generally 
at  the  low  free- stream  level,  implying  that  the  high  compression  ratio  across 
the  shock  is  not  achieved  except  in  the  vicinity  of  the  interface  of  the  two 
layers . 

Because  of  the  thin-layer  approximations,  the  thickness-curvature 
effects  which  generally  appear  in  the  higher-order  boundary  layer  or  Navier- 
Stokes  equations  do  not  appear.  Also,  the  slip  and  temperature -jump  effects 
can  be  separately  treated  as  higher-order  effects  of  the  shock-layer  theory. 
The  system  of  partial  differential  equations  governing  the  shock  layer  is  of 
parabolic  type,  for  which  a  forward  integration  in  the  downstream  direction 
is  possible.  The  corresponding  system  for  the  shock-transition  zone  is  re¬ 
ducible  to  one  of  ordinary  differential  equations,  very  similar  to  that  of  the 
one -dimensional  shock  wave.  This  system  gives  rise  to  a  set  of  modified 
shock  conservation  relations  which  account  for  the  transport  efforts  immedi¬ 
ately  behind  the  shock.  According  to  these  modifications,  the  total  enthalpy 
and  tangential  velocity  will  change  across  the  shock  in  a  manner  similar  to 
the  surface  slip  and  temperature  jump  across  the  Knudsen  layer  near  the 
body  surface.  With  these  shock  relations,  which  do  not  involve  any  detail 
of  the  shock-wave  structure,  the  flow  field  in  the  shock  layer  can  be  deter¬ 
mined  independently  of  the  shock-transition  zone.  An  interesting  general 
property  of  the  present  formulation  is  that  the  solution  always  has  the  correct 
values  for  skin  friction  and  heat-transf er  rate  in  the  f r ee-molecule  limit 
(for  unit  accommodation  coefficients). 
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The  shock-wave  structure  can  be  determined  by  simply  integrating  the 
ordinary  differential  equations  with  boundary  conditions  at  the  inner  edge 
furnished  by  the  shock-layer  solution.  The  analysis  provides,  in  essence, 
a  generalization  to  include  downstream  nonuniformity  under  the  condition 
5 s/a  44  I  . 

Unlike  the  existing  Newtonian  shock-layer  theories  for  inviscid  flow 

around  blunt  bodies,  the  present  analysis,  which  treats  the  tangential  pressure 

gradient  as  a  principal  term,  possesses  reasonably  good  accuracy,  as  has  been 

revealed  in  the  comparison  with  Van  Dyke's  inviscid  solution^  for  a  paraboloidal 

shock.  The  rather  detailed  agreement  with  the  merged-layer  solution  of 

40 

Levinsky  and  Yoshihara  confirms  the  validity  of  the  two-layer  concept  under¬ 
lying  the  present  formulation,  it  also  shows  the  surprising  accuracy  of  the  thin- 
layer  approximation  even  for  a  monatomic  gas. 

The  method  for  solving  the  dissimilar  problem  of  the  shock  layer  has 
been  developed  mainly  for  the  viscous  layer  and  the  merged-layer  regimes. 

It  provides  step-wise  integration  in  the  downstream  direction,  using  a  finite- 
difference  approximation.  The  difference  equations  are  singular  at  the  axis 
of  symmetry  and  at  the  body  surface.  On  account  of  the  singularity,  it  is 
necessary  to  develop  a  series  solution  in  the  vicinity  of  axis  of  symmetry 
to  produce  valid  initial  data  at  a  station  downstream  of  the  axis  of  symmetry. 
Because  of  the  singularity,  it  is  also  necessary  to  resort  to  the  use  of  an 
implicit  scheme  for  the  difference  approximation,  which  is  free  from  the 
usual  handicap  of  the  stability  and  convergence  requirements.  With  this 
implicit  scheme,  and  allowing  the  truncation  error  of  an  order  of  one  per¬ 
cent,  solutions  to  the  shock-layer  problems  can  be  brought  within  the  normal 
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capacity  of  a  standard  high-speed  digital  computer  (IBM  704). 

With  the  method  of  analysis  developed,  the  present  theory  is  applied 
to  analyze  flows  over  a  wedge,  cone,  paraboloid  and  hyperboloid  for  a  per¬ 
fect  gas  of  constant  specific  heat.  Generally,  the  solutions  obtained  reveal 
substantial  changes  in  total  enthalpy  and  tangential  velocity  across  the  shock 
in  the  low  Reynolds  number  regime  corresponding  to  £  1C  ^  0(1).  They 
provide  a  smooth  transition  from  the  boundary-layer  theory  to  free -molecule 
limit  (for  unit  accommodation  coefficients).  The  velocity  and  temperature 
gradients  near  the  body  surface  increase  with  distance  from  the  stagnation 
point  or  apex,  indicative  of  boundary-layer  developments. 

For  the  nonslender  wedge  and  cone  studied,  the  local  heat-transfer 
rate  and  skin  friction  are  found  to  be  governed  by  a  single  variable 
16  ~  ’  which  contains  the  gas  rarefaction  as  well  as  surface 

incidence  effects.  A  remarkable  feature  revealed  by  the  study  is  that  the 
velocity  profile  can  be  represented  very  closely  by  the  classical  boundary- 
layer  solution  (at  ib  3),  even  though  the  velocity  immediately  behind  the 
shock  is  quite  far  from  the  inviscid  value  and  the  viscous  layer  is  still  a 
major  fraction  of  the  shock  layer. 

For  the  problems  of  paraboloidal  and  hyperboloidal  shapes,  five 
successive  values  of  «2  are  considered  in  the  range  of  =0.3-  30, 

encompassing  the  viscous  layer  and  the  merged-layer  regimes.  For  these 
examples,  Tw/T0  =  1/8,  Pr  =  0.71  and  £  =  1/8.  Except  for  the  case  of 
lowest  Reynolds  number,  the  distributions  of  the  thickness  of  the  shock  layer 
around  the  body  are  relatively  insensitive  to  the  change  of  Reynolds  number. 
Comparison  with  the  inviscid  solution  shows  however  that  the  viscous  cases 
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studied  have  a  thinner  shock  layer  near  the  nose  region  due  to  strong  cooling 
by  heat  conduction  to  the  low-temperature  surface. 

The  results  obtained  reveal  interesting  heat-transfer  characteristics 
of  a  smooth  blunt  body  at  low  Reynolds  numbers.  The  ratio  of  local  heat- 
transfer  rate  at  a  fixed  station  <£  to  the  maximum  rate  at  the  nose  decreases 

from  the  free -molecule  distribution  with  increasing  density.  A  reversal  of 

2, 

this  trend  occurs  however  at  K  ca  10,  corresponding  to  the  viscous-layer 

regime.  This  reversal  can  be  explained  by  the  vorticity  effect.  One  may 

recall  from  the  discussion  that  the  variations  in  the  heat-transfer  ratio  in  the 

range  of  «  =  3-30  are  remarkably  small.  This,  together  with  the  fact  that 

the  vorticity  interaction  is  itself  not  a  large  effect,  suggests  that  there  is  a 

very  wide  range  of  Reynolds  number  (corresponding  to  the  boundary-layer 

through  the  incipient-merged-layer  regimes)  in  which  the.  heat-transfer 

distribution  on  a  smooth  blunt  body  remains  essentially  unchanged.  Similar 

behavior  has  also  been  observed  for  the  skin  friction  characteristics. 

In  the  study  of  nonequilibrium  dissociation  and  species  diffusion  in 

the  stagnation  region,  an  example  is  carried  out  for  an  idealized  binary 

mixture  of  "air  molecules"  and  "air  atoms",  and  a  cold  body  surface  at 

3 

near  -  satellite  speed,  comparable  to  the  previous  analysis  of  Chung.  The 
results  reveal  that  the  "concentration  jump"  across  the  shock  resulting  from 
the  upstream  diffusion  of  the  atomic  species  is  appreciable  even  in  the 
viscous -layer  regime  where  changes  in  the  total  enthalpy  and  tangential 
velocity  across  the  shock  may  be  neglected.  In  the  merged-layer  regimes, 
as  a  result  of  the  combined  effect  of  surface  cooling  and  enthalpy  drop  across 
the  shock,  the  gas-phase  chemistry  is  practically  frozen.  Even  in  the 
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viscous-layer  regime,  dissociation  takes  place  only  in  the  outer  part  of 
the  shock  layer  near  interface.  In  these  low-density  regimes,  the  surface 
catalycity  has  a  predominant  influence  over  the  composition  in  the  inner 
portion  of  the  shock  layer,  as  is  expected.  However,  there  is  very  little 
change  in  the  temperature  profile  associated  with  change  in  surface  cata¬ 
lycity,  because  of  the  frozen  gas-phase  chemistry  in  the  inner  part  of  the 
shock  layer. 

The  discussion  has  also  brought  out  the  importance  of  pressure 
diffusion  in  shock  layers  removed  from  the  stagnation  region. 
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APPENDIX 


In  the  following,  important  details  of  the  present  method  of  numerical 
solution  to  the  shock-layer  equations  described  in  Section  6  will  be  provided. 


A.  1  The  System  of  Ordinary  Differential  Equations  Governing  the 
Coefficients  U,  ,  ©,  ,  and  p'-x  . 


For  blunt  bodies  with  finite  nose  radii,  Z 

k-^4- 

4  A 


+ 


is  expressible  as 


The  function  U  ,  <9  and  -p-  may  assume  the  form 

U  -  U0+  \u,  + . 

0  =  0O  +  A  &!  +  •  •  ’ 

-py  ~  i  +  A  -fx,  +  \  p'^  4-  ....  . 


The  system  of  equations  governing  u.o  ,  0o  and  -jE-  have  been  given  in 
Section  6.  2  of  Part  II  The  system  governing  LLt  ,  0,  and  are 

provided  below. 
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u,(o)  =  @,(o)=o 
|  frO)  =  o 

|«/0  =  uj-u0)(i>~ k)-(^-£L  -|1  a; 

U  c)  -  o  -  - 1  ♦  -y  -  5- ('-  “J 

A.  2  The  Implicit-Difference  Formulation  of  the  Shock-Layer  Equations 

The  complete,  final  form  of  the  finite-difference  formulation  of  the 
shock-layer  equations,  based  on  the  implicit  scheme,  are  given  below.  For 
each  point  e  ,  there  are  two  algebraic  equations 

A  u-ec  +  B  +  C  -  Ox 
D  +  E  @e  +  F  ©f  =  J 

which  relate  the  values  of  u  and  @  at  the  point  e  to  the  values  at  the 
neighboring  points  cC  and  -f  . 

The  subscripts  refer  to  grid  points  specified  in  the  sketches  below  and 
in  Fig.  6.1. 


(A.  1) 


(b) 


$  = 1-0 


The  coefficients  A,  B,  C,  D,  E,  F,  C,  and  J  in  Eq.  (A.l)  are: 
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between  f  and  ^  =  1  by  application  of  Eq.  (A.l)  is  to  be  solved 

with  the  boundary  conditions: 
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Figure  7.9  DISTRIBUTION  OF  SURFACE  HEAT-TRANSFER  RATE  AROUND 
A  PARABOLOID  AND  HYPERBOLOID 
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line,  based  on  the  present  method,  after  correction  is  in  g*j®d 
agreement  with  th#  w1'  M  \ '.n  ihvke  and  ?w.-tga rt 

p  1 05  Exchange'  - - - — y*  and  - r — --  i  j  Fig  i 

i-t„***F  Jao'Un/ 

P  109  Figure  7  ,i  a  extend  the  arrow  fr  m  “3  0"  down  to  the  third 
curve  r.  r'-onted  -  t  m  p  il'.e  r  pit.  »- r  right 

p  J19  Figure  8  2  c .  the  •  »■.«*•  ( mm  "K"  '88"  sh<o«jd.-i  be  pointed 

in  the  third  ,  rve  ;  »  m  the  b  tt  •  > ? a . 


